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Abstract 



In this thesis we examine three different models in the MSSM context, all of which 
have significant supergravity anomaly contributions to their soft masses. These mod- 
els are the so-called Minimal, Hypercharged, and Gaugino Anomaly Mediated Super- 
symmetry Breaking models. We explore some of the string theoretical motivations 
for these models and proceed by understanding how they would appear at the Large 
Hadron Collider (LHC). Our major results include calculating the LHC reach for each 
model's parameter space and prescribing a method for distinguishing the models after 
the collection of 100 fb" 1 at y/s = 14 TeV. AMSB models are notorious for predicting 
too low a dark matter relic density. To counter this argument we explore several 
proposed mechanisms for non-thermal dark matter production that act to augment 
abundances from the usual thermal calculations. Interestingly, we find that future 
direct detection dark matter experiments potentially have a much better reach than 
the LHC for these models. 
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Introduction 



Particle physics is at an exciting point at the time of writing this thesis. The 
long-awaited start of the Large Hadron Collider (LHC) era has finally begun with the 
first collection of data at the world-record breaking collision energy of 7 TeV c.o.m. 
in March of this year. After the first two years the experiment is planned to run 
at the full design luminosity accumulating ultimately 100-1000 fb" 1 at a tremendous 
y/s = 14 TeV. It is expected, or at least hoped, that the LHC will shed light on 
important mysteries of the Standard Model (SM) of particle physics by allowing us 
to detect particle states that have yet to be observed at the Tevatron. The most likely 
of these is the Higgs boson, the remaining piece of the SM, which is thought to be 
responsible for the spontaneous breaking of the electroweak symmetry. The presence 
of the Higgs boson in the SM is itself a strong theoretical motivator that other heavy 
particle states exist that can also be discovered at the LHC. One very well-motivated 
class are the supersymmetric (SUSY ) particle states whose masses lie in the TeV 
range. The discovery of these particles at the LHC would have deep implications for 
the nature of space-time. Having the potential to discover new physical states such 
as these puts us at a truly unique and exciting time. 

On the other hand particle physics is also merging increasingly with cosmology. 
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The energy density of the Universe is known to be comprised of Dark Energy (71%), 
normal baryonic matter (4%), and a non-luminous form of matter known as Dark 
Matter (DM) which comprises roughly 25% of the energy density. If DM is consid- 
ered to be comprised of particles, it must be massive to account for the relic abundance 
and it must be cold enough to allow for structure formation on large scales. For these 
reasons, there are no good candidates for DM in the SM. There is, however, a partic- 
ularly good DM candidate in supersymmetric theories in which R-parity is conserved. 
The particle is the lightest neutralino, and the search for it is an important priority 
at the LHC and in DM experiments around the world. 

The outline of this thesis is as follows. The remainder of the Introduction is in- 
tended to provide background for subsequent chapters. We will first give a very brief 
introduction to the Standard Model (SM) and simultaneously attempt to motivate 
supersymmetry. Next, we will introduce the idea of the scalar superfield before mov- 
ing on to describe the Minimal Supersymmetric Standard Model (MSSM)Q Once the 
fields and interactions of the theory are established, we turn our attention to how the 
scalar superpartners of the SM acquire mass at the weak scale. These quantities are 
usually given by a combination of low-energy constraints (e.g, electroweak symmetry 
breaking parameters, etc.) and by high-scale physics (M GUT , M string , etc.) that re- 
lates to physics at the low scale (M wea k) through renormalization group effects. For 
this thesis, the latter quantities are, in part, due to anomalies that are present in 
supergravity theories, and a rather technical derivation of these will be given Section 
1.4| Theories in which contributions of this type are important will be referred to as 
"Anomaly- Mediated Supersymmetry Breaking (AMSB) models". The Introduction 
concludes with descriptions of the computational tools used in this work and should 
be used for reference in the chapters that follow. 

In Chapters [2] - [4] we will examine three different models of the AMSB type. The 

1 Here, the SM and MSSM are best viewed as the low-energy effective theories of a high-energy 
framework, i.e., supergravity. 
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first is the prototypical model known as the minimal AMSB model, and is "min- 
imal" in the sense that it has the minimum ingredients to be phenomenologically 
viable. The second model is referred to as "Hypercharged AMSB" (HCAMSB) be- 
cause it pairs AMSB mass contributions with [/(l)-contributions at the string scale. 
The third model is actually a class of string theory models with specific high-scale 
boundary conditions and a rather low string scale, referred to collectively as "Gaug- 
ino AMSB", or "inoAMSB" for short. Each of these chapters gives a full theoretical 
explanation of the model being considered. They also give a full analysis for LHC 
physics, that is, we describe all renormalization group effects, compute weak-scale 
parameters, run signal and SM-background simulations for pp collisions at y/s = 14 
TeV, and calculate the reach of the models' parameters for 100 ftT 1 of accumulated 
data (one year at full design luminosity). 

It would be greatly insufficient to focus only on collider physics searches since 
future cosmological data will be precise enough to be competitive. In Chapter [5] we 
consider the question of dark matter (DM) in AMSB models where we take the light- 
est neutralino to be our DM candidate. AMSB models notoriously yield too little 
thermally-produced relic abundance to account for the measured DM. However, it is 
expected that other heavy particles present in the early universe could increase the 
DM abundance in AMSB models non-thermally . In this chapter we present several 
such scenarios and assume the total DM abundance can be accommodated. After 
a review of DM theory and direct /indirect detection experiments, we calculate im- 
portant rates for AMSB DM cosmology for all of the models in Chapters [2] - |4j We 
supplement all rates with detailed physical explanations. We also find the reach for 
each model and compare the results with LHC expectations. 

Finally, we conclude with a general overview of the results in Chapter [6] 
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1.1 Standard Model and Supersymmetry 



1.1 Standard Model and Supersymmetry 

So far most experimental data in HEP experiments is described by the Standard 
Model of particle physics and this has been the case for 30+ years^J The basic facts 
of this model are given in this section with the goal in mind of extending it to include 
supersymmetry. 

The Standard Model is a collection of three identical generations of spin-^ fermion 
fields, spin-1 vector bosons, and spin-0 scalar Higgs fields. The fermions interact 
through the exchange of spin-1 vector bosons that arise through the gauge invariance 
of both abelian and non-abelian interactions. The interactions of the particles can be 
trivial or non-trivial under each of the SU(3) C , SU(2)l, and U(1)y rotations, where 
"c" stands for the color-charge and "L" for the weak-isospin, and "Y" is the hyper- 
charge of the interacting particles. 

Those particles that interact through SU(3) C interactions are the spin-| fermions 
known as quarks and the interaction is mediated by spin-1, color-charged octet of 
gluon fields. Fermions that are unaffected by SU(3) C interactions are known as lep- 
tons. Left-handed fermions interact in very specific pairs (or doublets) and the me- 
diators of the interactions are the spin-1 triplet of W^-bosons. These interactions do 
not apply to right-handed fermions, and they are the only sources of flavor-change 
in the SM [8] as is described, for example, by Cabibbo-Kobayashi-Maskawa mixing 
for quarks. And finally, all fermions are charged under abelian U{1)y rotations, and 
these interactions are mediated by the vector boson. 

The inclusion of a Higgs particle is necessary to give mass to both the weak sector 

bosons and matter fields. The Higgs ($ = (^o)) is thought to be a dynamic field with 

important exceptions not described by the SM are ^-oscillations, dark matter, dark energy, and 
the baryon asymmetry. 
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quartic self-interactions, thus its contribution to the Lagrangian is 

Cm g9 s = \d^\ 2 + ^ - A($ t $) 2 . (1.1) 

By allowing $ to be a complex doublet that transforms non-trivially under the elec- 
troweak symmetry, SU(2)l x U(l), the covariantized action mixes the Higgs with 
weak and hypercharge bosons [90] . When the neutral Higgs acquires a VEV, the 
SU(2) x U(1)y symmetry is still a good symmetry of the Lagrangian, yet its generators 
no longer annihilate the vacuum. The symmetry is said to be broken spontaneously , 
and the remnant symmetry is the U(l) em of electromagnetism. In the symmetry 
breaking process, massless degrees of freedom are produced which are subsequently 
transformed into the longitudinal modes of the electroweak bosons (rendering them 
massive). Also upon acquiring a VEV, the Higgs gives mass to the fermion fields 
through its Yukawa couplings to matter. 



The fundamental particles of the Standard Model are summarized in Table 1.1 
along with the matter and Higgs quantum numbers. The table tells how the par- 
ticles transform under a given symmetry. For example, left-handed fermions form 
doublets transforming non-trivially under SU(2), right-handed fermions are singlets 
under SU(2), and quarks transform as a triplet of color charge under SU(3). Not 
shown in the table are two other generations of fermions, identical to the one shown 
with exactly the same quantum numbers, but with larger masses. Also, anti-particles 
have not been included. 

The SM is completely consistent but there are several reasons why it cannot be 
the complete description of nature. A list of some, but not all, of these reasons are 
given here: 

• hierarchy problem - radiative corrections to the scalar (Higgs) mass terms are 
quadratically divergent due to gauge and fermion loops, but unitarity arguments 
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Fermion Flavor 


Symmetry/ Quantum #s 
SU(3) C SU(2) L U(1) Y 


leptons \ e i l 

e R 

quarks KdJ L 
ur 


12 -1 
11 -2 

3 2 | 
3 1 | 
3 1 | 


boson v * o; 

(1,2,1) 


r gluon weak hypercharge 
gauge I W* B, 

L A=l---8 i = Q,+,- 



Table 1.1: Matter, vector boson, and Higgs fields in the Standard Model. 



require the mass to be constrained to less than a few hundred GeV [20]. If the 
SM is to be considered an effective field theory below a high-scale A, then the 
Higgs mass could be subject to excessive and unnatural fine-tuning. 

• Dark Matter (DM) is not yet included in the SM. 

• Gravitational interactions are not present in the SM. 

It is ideal to have a single framework that can address these issues, and in this the- 
sis supersymmetry is the adopted solution. A supersymmetry transformation acts 
on bosonic state to form fermionic states and vice versa. The hierarchy problem is 
famously solved in supersymmetric theories because the fermion loop corrections to 
the Higgs mass are accompanied by bosonic corrections. The new loops are also 
quadratically divergent but generally appear with opposite sign. In the case of 
unbroken supersymmetry, where fermions and bosons have equal mass, the fermionic 
and bosonic contributions precisely cancel one another to all orders of perturbation 
theory [20j[82j[101j. Even in the case of broken supersymmetry, the divergent con- 
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tributions to the Higgs mass are at most logarithmic and, therefore, not severely 
destabilizing. 

Supersymmetry also provides a framework for addressing the remaining two issues 
in the list above. In particular, we will consider local supersymmetry at a very high 
scale that naturally incorporates gravity. It will be necessary for the supersymmetry 
to be (spontaneously) broken to agree with phenomenology. When we define the 
minimal supersymmetric model and renormalize the model parameters at the weak 
scale, we will encounter natural EW symmetry breaking, non-quadratically divergent 
scalar masses, and supersymmetry provides several DM candidates. 

Before developing the superfield formalism in the next section, it is important 
to discuss some technical details regarding why, if supersymmetry is to explain the 
short-comings of the SM, it must be a broken symmetry. Since the supersymmetry 
generators transform bosons into fermions, they must be fermionic and therefore obey 
anti-commutation relations. When these relations are combined with the Poincare 
algebra, the closed (graded) algebra that results is known as the super-Poincare al- 
gebra. The squared-momentum generator, (P M ) 2 , is a Casimir of the super-Poincare 
algebra [81], and so supersymmetric fermion-boson pairs are expected to be degen- 
erate. If this was a rigid requirement, supersymmetry would already by excluded by 
the fact that no partners of the SM particles have been observed with identical mass. 
Supersymmetry necessarily has to be a broken symmetry to be phenomenologically 
viable. 

Viable supersymmetric models must assume breaking in a way that does not re- 
introduce quadratic divergences. In order for this to occur it is necessary for the di- 
mensionless couplings of the theory to be unmodified by supersymmetry breaking, and 
that only couplings with positive mass dimension are included in the supersymmetry 
breaking potential. Supersymmetry that is broken in this way is said to be broken 
softly. This reduces the number of additional terms that can be included in the soft 
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supersymmetry-breaking Lagrangian, C so ft, which are: linear (gauge singlet), bilinear 
(including masses), and trilinear (A-term) scalar interactions, and bilinear gaugino 
mass terms. 



In this section we briefly review some aspects of the superfield formalism. We will see 
how boson-fermion pairs are embedded into super-multiplets. Also in this section, 
key ingredients to supersymmetric theories are defined, and the notation to be used 
in describing the minimal supersymmetric standard model of the next section are 
established. 

Superfields place boson fields in the same multiplet as their fermionic partners. It 
is convenient to introduce anti- commuting Grassmann variables arranged as Majorana 
spinors, 9 (6), that multiply the fields in order to put them on the same footing. Fields 
in a super-multiplet that are not multiplied by Grassmann variables are referred 
to as the "lowest" component, those with one Grassmann variable are the second 
component, etcj^j It is the lowest component that determines the type of superfield, 
i.e., scalar, spinor, vectoiQ, etc. 

A general scalar superfield has many components as seen in 



1.2 Superfields 





3 For more discussion on Grassmann variables, see [55] |102) . 

4 We only consider scalar superfields here. For more information, see |20) . 



s 



1.2 Superfields 



but this representation is not irreducible. Fortunately supersymmetry allows for chiral 
representations of superfields. For example, there is a representation where the scalar 
component and the left-chiral spinor transform into one another without mixing with 
the corresponding right-handed fields. A left chiral scalar superfield is of the form 

S L (x, 9) = S{x) + iV29ip L {x) + m L F(x) (1.3) 

where the lowest (Grassmann) component is a scalar, the second component is the 
partner fermion, and F is an auxiliary field required to balance off-shell degrees of 
freedom^} Similarly a right chiral scalar superfield is of the form 

S R ( X , 6) = S(tf) - iV26ip R (x j ) + i66 R F(x j ), (1.4) 

but will frequently be recast as the conjugate of a left-chiral scalar field. In short- 
hand, chiral scalar superfields can be referred to by their components, (S,ipL,F). 
To be clear, the scalar component is not a spinor and does not have helicity. It 
contains the annihilation operator of the superpartner of a chiral fermion. Under 
supersymmetry transformations the components of the left-chiral superfield transform 
into one another: 

5S = -iV2atfj L (1.5) 
6ip L = -y/2Fa L + y/2$Sa R (1.6) 
5F = iV2a0ifj L . (1.7) 

It is intriguing that the F-term of the superfield transforms as a total derivative 



(Equation (1.7)) because field combinations that transform as total derivatives bring 



Additionally x^ = x^ + idj^j^O, but we will not focus on technical details of this sort. 
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dynamics to the action. This property of the F- term is true even of products of chi- 
ral superfields because, as can be shown, they are themselves chiral superfields. The 
general polynomial of left-chiral superfields is another left-chiral superfield known as 
the superpotential , and its F-term (00 l- component) must be appended to the La- 
grangian. In a renormalizable theory, the superpotential is at most a cubic polynomial 
by dimensional analysis. 

It can also be shown that the (#75#) 2 -component, or "D-term", of a general su- 



perfield (Equation (1.2)) transforms as a total derivative under supersymmetry trans- 
formations . The product of chiral superfields, however, does not have an interesting 
D-term because it is already a total derivative, as in 



Sl 3 l(0j 5 0) 2 V 2 S. (l.i 



This term is automatically zero in the action and does not produce any dynamics, 
and is for this reason that the superpotential does not contribute D-terms. However, 
a polynomial of mixed chirality can give important D-term contributions to the La- 
grangian, and this function is referred to as the Kahler potential. It is important to 
note that the Kahler potential is at most quadratic by renormalizability and real by 
hermiticity of the Lagrangian [20J. It is then taken to have the form 

TV 

K=J2 4*4%, (1.9) 

the D-term of which is also included in the Lagrangian as with the F-term contribu- 
tions. 
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1.3 Minimal Supersymmetric Standard Model 

We are now ready to build the minimal extension of the Standard Model (SM) that 
incorporates supersymmetry, the Minimal Supersymmetric Standard Model (MSSM). 
That is, we seek the minimal extension of the SM that includes broken supersymmetry, 
and that is both phenomenologically and theoretically safe. In order to do this, it 
is first assumed that the theory will have the SM gauge symmetry group: SU(3) C x 
SU(2)l x U(l). Thereafter, all gauge and matter fields of the SM must be promoted 
to superfields. The matter superfields must be L-chiral fields as required by the 
superpotential. Thus, for each chirality of every SM fields, we will choose there will 
be one Ze/t-chiral superfield assigned. 

Extending the SM Higgs field to a superfield will add a partner fermion, a higgsino 
(h). Having only one extra fermion re-introduces U(l)y and U{1)ySU{2) 2 l gauge 
anomalies [HS] that are canceled successfully in the SM. If instead there are two 
Higgs doublets in the theory, with opposite hypercharges, Y = 1 and Y — — 1, 
their fermionic partners have the correct quantum numbers to satisfy the anomaly 
cancellation. Furthermore, a single Higgs doublet is not allowed in the MSSM because 
the lower- component fermions of the SM weak doublets would receive their mass from 
the conjugate of the Higgs, which is a right-chiral superfield. Interactions of right- 
chiral superfields are forbidden in the superpotential [20] , and we are forced to accept 
at least two Higgs doublets into the MSSM. We denote by H u the Higgs doublet of 
superfields that is associated with the mass of Y = 1 fermions, and Hd associated 
with mass of Y = — 1 fermions. The matter and Higgs superfields of the MSSM are 
shown in Table 11.21 

Next it is necessary to define the superpotential of the theory. The superpoten- 
tial, denoted by W, contains SU(2) x U(l) invariant combinations of chiral matter 
superfields and Higgs fields. The matter fields are coupled to Higgs fields through 
Yukawa coupling matrices, while the /i-term couples H u and H d . In the MSSM it is: 
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Field 


SU(3) C SU(2) L U{1) Y 


E c 

Q=(t) 

u c 

D c 

TT (ht\ 

&<=(%) 


12-1 
112 

3 2 § 
3* 1 -f 
3* 1 | 

1 2 1 

1 2* -1 



Table 1.2: Matter and Higgs superflelds in the MSSM. 

W = fiH:H da + [$«)MHjj C j + (ZdaQtHdafrj + {f e )i j L a i H da E*\, (1.10) 
»,i=i,3 

where a is an SU{2) index and c stands for conjugation. 

This superpotential is not completely general because, in its construction, terms 
that would lead to baryon (B) and/or lepton (L) number violation have been carefully 
omitted. There are renormalizable terms that could be added that are gauge and 
supersymmetrically invariant, but the presence of these terms would have physical 
consequences that are highly constrained by experiment (for instance, proton decay). 
The omission of these operators is made possible by imposing a discrete symmetry 
on the superpotential, known as R-parity. When 

R = (-lfP-L)+2. (l.H) 

is conserved (s is the spin of the state), renormalizable B- and L- violating interactions 
will be forbidden . Furthermore, conservation of R-parity leads to three important 
consequences [99] : 
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1. superpartners will always be produced in pairs at colliders; 

2. the decays of the SM superpartners (including extra higgsinos) produce an odd 
number of the final state lightest SUSY particle (LSP), which for our purposes 
this is a neutralino; 

3. the LSP is absolutely stable and therefore may be a good dark matter candidate. 

We accept R-parity as part of the definition of the MSSM and should stay mindful 
of these consequences in the coming chapters. 

The final step in constructing the MSSM is to include all gauge-invariant soft 
supersymmetry-breaking terms into the Lagrangian. These terms are thought to 
arise from the interactions between MSSM fields and a "hidden" sector where SUSY is 



broken (see Section 1.4). These terms raise the masses of supersymmetric partners of 
the SM fields and are needed for electroweak symmetry breaking (see next subsection). 
The MSSM soft Lagrangian is [20] 



L\ m lijLj + 4i m E;/iy + m H u \ H u\ 2 + m 2 Hd \H d \ 2 ] 



- l - [M x AoAo + M 2 X A X A + Mfg B ~g B ] 

- - [M^AoTsAo + M' 2 \ Als \ A + M^ b 159b\ 

+ [(c u ) tJ e ab Q«H* d b rf Rl + (c d ) tJ Q«H* u J Rj + (c e ) tj L?H* u J R] + h.c] 

+ [bH a u H da + h.c]. (1.12) 



In short, the parameters above are the scalar mass matrices of lines 1 & 2, the gaugino 
mass terms of lines 3 & 4, the trilinear A-term couplings from supersymmetry breaking 
in line 5, trilinear terms of line 6 & 7. It is seen that these terms give the MSSM 
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an extremely large number of parameters. We will see that in AMSB-type models 
the number of parameters will be reduced from the 178 in the soft Lagrangian above 
to just a few. This is one of the many attractive features of the class of models 
considered in this thesis. 



Electroweak Symmetry Breaking in the MSSM 

As in the Standard Model, the electroweak symmetry is spontaneously broken by 
minimizing the potential in the scalar sector. We have seen that due to supersymme- 
try the scalar sector includes much more than just the Higgs particle. The potential 
is extended now to include effects of all matter scalars along with all possible effects 
that originate in SUSY breaking. The scalar potential is then the sum of the various 
terms: 

Vmssm = V f + V d + V soft (1.13) 

with 



dW 



s=s 



and 



;i.i4) 
;i.i5) 



and V so ft comes from lines 1, 2, 5, 6, and 7 of Equation (1.12). As in the SM, the 



Higgs VEVs should be responsible for the breaking. In the usual manner, gauge sym- 
metry allows the VEV of H u to be rotated to its lower neutral component. Upon 
minimizing with respect to the other component of H u , it is necessarily so that {h~[) 
= [20] . Then only the potential of the neutral Higgs fields needs to minimized in the 
breaking of electroweak symmetry. In this way, provided no other scalars are allowed 
to develop VEVs, only charge- conserving vacuua can occur in the MSSM. 

It will be important to understand how the electroweak symmetry is broken prop- 



14 



1.3 Minimal Supersymmetric Standard Model 

erly in order to put constraints on parameters in a model. Schematically the mini- 
mization requires the following: 

1. the potential is extremized with respect to both h° u and h° d and their conjugates 
through its first derivatives, i.e. 



2. to ensure that EW-breaking occurs the origin must be a maximum, i.e., the 
determinant of the Hessian should be negative there, and this imposes the con- 
dition 



3. and finally, in order that the potential is bounded from below the (D-term) 
quartic terms must be non-zero (positive at infinity) and this results in the 



Electroweak symmetry is broken properly when these conditions are met because a 
well-defined minimum develops that does not break charge. It is customary to define 
the ratio of the Higgs VEVs as 

tan/?=(^), 

and to recast the important potential minimization conditions as 



8V 
dh? 



(and conjugates) = 0; 



(B^ 2 >(m 2 Hu +^ 2 )(m 2 + // 2 ); 



condition 



2\Bfi\ <m 2 Hu +m 2 Hi + 2^. 



(m 2 Hu + m 2 Hd + 2/i 2 )sin2/3 



and 



(1.16) 




2 

m H d -m 2 Hu (Un f3) 2 M 2 



(1.17) 



tan 2 /3 - 1 2 
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where Mz = 9 t, g (v* + t>J) is the tree-level result for the Z° mass when the effects 
of EW-breaking on gauge bosons are considered. 

In the SM, the required shape of the potential is achieved through the inclusion of 
a tachyonic scalar that transforms non-trivially under the EW symmetry and acquires 
a VEV. The same is true in the MSSM, with the exception that the potential is gen- 
erated naturally at the weak scale through RGE effects rather than put in by hand. 
As SUSY parameters are evolved from the GUT scale down to the weak scale, the 
large top Yukawa coupling drives the up-type squared Higgs mass to negative values, 
thus triggering what is know as radiative electroweak symmetry breaking (REWSB). 



The fj, value calculated through Equation 1.17 is purely phenomenological. In 
actuality, it is difficult to understand why \i should be so small considering it appears 
in a supersymmetric superpotential term and should naturally be of order the (high) 
supersymmetry breaking scale (perhaps of order the Planck scale, M P i) [38], but such 
high values would destroy the mechanism of electroweak symmetry breaking. This is 
known as the "/i-problem" . Typically the low \i value is assumed to arise by some 
other mechanism, most commonly that of Guidice and Masiero [59J. 

The procedure described above is true in general for minimizing the scalar po- 
tential, however there are important cases where radiative corrections need to be 
considered. For example, without radiative corrections, the MSSM would already be 
excluded by Higgs mass bounds. For the sake of brevity we do not dwell on these 
corrections, but their importance is noted. 

Now that we have seen the MSSM and its parameters, we can determine the physi- 
cal mass eigenstates that are important for any type of phenomenology. The following 
subsections give a general overview of the contributions to the mass matrices at the 
weak scale in the various sectors. 
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Neutralinos and Charginos 

The spontaneous electroweak symmetry breaking leads to mixing of fields with the 
same electric charge, spin, and color quantum numbers. The spin-|, color- neutral 
fermions mix and come as wino-bino-higgsino mass eigenstate combinations. Those 
combinations that are neutral are called neutralinos and those that are charged are 
charginos. 

The Lagrangian mass terms for the neutral fields can be written as 



C 



1 



neutralino 



neutralino 



'1.18) 



where \1/ T = hj, W 3 , b) contains the neutral up- and down-type higgsinos, wino, 
and bino respectively. The tree-level mass matrix^ for this sector is 



/ 



M, 



neutralino 







/' 

gy-u, 

\ V2 





g'v<i 

V2 



gvu 
gvd 





g vu 

V2 
V2 



\ 





M l J 



1.19) 



After diagonalizing the mass matrix we find the mass eigenstates, Zi, in linear com- 
binations of the higgsinos and gauginos, i.e., 



Vi(i)h° u + v i{2 )h° d + v i(3 )W 3 + Dj( 4) 6 



1.20) 



Note that the neutralino mass matrix is hermitian and will have real eigenvalues. 
Charginos are linear combinations of the charged gauginos, W , and charged 



6 The actual mass matrix has 1-loop contributions, the most important of which on the upper-left 
diagonal. 
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higgsinos, h u and h d . They appear in the Lagrangian as 



chargino 



MMchargedPL + M^ harged P R )^ , 



1.211 



where the negatively charged, two-component field is $ T = ( Al +' Az ; P L h d — Pr/i+), 
and the tree-level matrix is 



M 



M 2 -gv d 



charged 



\ 



-gvv 



(1.22) 



7 



The diagonalization of M. charged is more involved than in the neutral case as it requires 
two unitary matrices, U and V, and this procedure can be found in [20J. In the end, 
there are two mass eigenstates for each charge in the combinations 



Wf = U i(1) W- + U m h~ 
w? = v i{1) w + + V i{2) h+. 



;i.23) 

1.24) 



The physics mass eigenvalues will depend on /i, M2,Mw, and tan j3. Note that the 
mass matrix is tree-level and loop corrections are not shown here. 

In general, there may be a high degree of mixing for any of the neutralino and 
chargino mass eigenstates. It turns out that of the parameters fi, Mi, and tan /3 
that contribute to the mixing, AMSB type models tend to have a very light M 2 value. 
This results in both the lightest neutralino and lightest chargino being wino-like, a 
fact that will have strong implications for LHC phenomenology (Chapters [2]- [4]) and 
Dark Matter (Chapter [5]) as we will see. 
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Sfermion Masses 

Unlike the SM, there are four mass contributions to the squarks and sleptons. They 
are: superpotential contributions that rely and Higgs VEV; generational soft masses 
from SUSY breaking; trilinear interactions with Higgs fields; and D-term contribu- 
tions. The result of these contributions is to mix left and right sfermions of the same 
flavor with 2x2 matrices that are straight-forwardly diagonalized. The procedure 
is not very illuminating for the current discussion, and the reader is referred to [20] 
[82"] for more discussion. In the end, the mass eigenvalues will depend on /i, tan (3, 
trilinear couplings, SM EW parameters, and the corresponding fermion mass. 

Gluino Mass 

The gluino mass is not tied to the EW symmetry breaking sector, and its presence is 
purely due to supersymmetry breaking and renormalization, and is parameterized by 
M3. The gluino must be a mass eigenstate because it is the only color octet fermion 
and SU(3)c is not a broken symmetry. In the next section we will see how this and 
other gaugino masses are generated by anomalous SUGRA effects. 



1.4 Supergravity Soft Terms 

The MSSM is considered to be an effective theory containing information about su- 



persymmetry breaking through the soft terms of Equation (1.12). In this section the 
origin of the soft terms is outlined. Most importantly, the explanation of how the 
supergravity anomaly imparts mass is explained. The class of models where these 
contributions dominate (or is comparable to) all other forms of soft mass generation 
are known as Anomaly Mediated Supersymmetry Breaking Models, and all of the 
models in this thesis are of this type. 
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It has already been remarked that the generators of supersymmetry transforma- 
tions are spinors, and it follows that the variational parameter of supersymmetry 
is also a spinor. When supersymmetry is "localized" and combined with the sym- 
metries of general relativity, the resulting theory is called supergravity [52j[102j[82j. 
Supergravity has a supermultiplet of gauge particles: the spin-2 graviton, and its 
super-partner, the spin-| gravitino, both of which are massless so long as supersym- 
metry is unbroken. However, when supersymmetry is broken spontaneously by the 
VEV of the auxiliary component of a "hidden" supermultiplet X, i.e., < Fx >, a 
massless goldstone fermion, or goldstino, is produced. This goldstino becomes the 
longitudinal mode of the gravitino thereby imparting mass, and this mechanism is 
called the super-Higgs mechanism due to its obvious parallels to SM Higgs mecha- 
nism. The gravitino mass, m 3 / 2 , is given by 



<Fx> n 

m3/2 = T^pT' (L25) 



and plays a key role in soft mass generation when supersymmetry breaking is com- 
municated to the MSSM particles. In the case of gravity-mediated supersymmetry 
breaking, MSSM fields interact with the hidden sector mainly through gravitational 
interactions. These interactions induce soft masses for the scalars that are suppressed 
by powers of Mpi and typically of the order 



m soft ~ < ff > - (1-26) 
Mpi 



When the MSSM is coupled to gravity, the effective soft Lagrangian below the Planck 



scale is given as in Equation ( 1.12 ) with coefficients that depend on powers of m so f t ~ 
m 3 / 2 . However these tree-level effects can be suppressed relative to the quantum 
anomaly contributions to be described in the next subsection, and we will see that 
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these contributions will require typically higher values of m 3 / 2 than is typical for 
gravity mediated scenarios. 



Anomaly Mediated Supersymmetry Breaking 

In this section we examine how soft masses are generated by anomalous supergrav- 
ity (local susy) effects^] Low-energy type II-B string theory enforces the following 
functional as the classical Wilsonian action for a generic super Yang-Mills coupled to 
gravity (see [55]|102] for more discussion): 

S = - 3 J d 8 zE exp{-^K($, $; Q, Qe 2V )} (1.27) 
+ ( J d 6 zS \W(Q, Q) + \f ab W aa W b a ] + h.c.) 
+ (higher order derivatives), 



where K is the Kahler potential, W the superpotential, f a b are the gauge coupling, 
and W a is the gauge field strength associated with the prepotential V. R is the chiral 
curvature superfield, E is the full superspace measure, and £ is the chiral superspace 
measure. 

The issue of the anomaly arises with the first term in the functional above because 
it is not in the standard form with Einsteinian gravity. It also has the improper form 
for kinetic metric, given by the second derivative of the Kahler potential. What is de- 
sired instead is to have the fields transformed out of the SUGRA frame (as it appears 



in (1.27)) into the more physical Einstein-Kahler frame that has Einsteinian gravity 
as well as canonical matter normalizations. The supergravity effects are disentangled 



7 The following description follows closely to the arguments of Kaplunovsky and Louis [6 8) and 
de Alwis 1461. 



21 



1.4 Supergravity Soft Terms 



via Weyl-rescaling of the metric, i.e., 

g-mn = Qmn " e~* K , (1.28) 

and by rescaling fermion fields by exponentials of K in order that fermions and bosons 
belonging to the same supermultiplet are normalized in the same way [102j[68j. 

The torsion and the chirality constraints of supergravity (SUGRA) are invariant 
under Weyl transformations, given by re- weightings with parameter e T and e r , where 
t is a chiral superfield. Some of these transformations are given below: 

E -> e 2 ( r+f )E, £^e 6r £, 
V^V, V a ^e T ~ 2f (V a + ...), 

Q^Q, & W a ^e~ 3T W a . 

The SUGRA action itself is not invariant under these transformations. However, we 
can make insertions of the "Weyl compensator", C, such that the action is invariant: 

r -v 2 - 1 

S=-3 d 6 z£ (—j- + 2R) CC exp{--K($, $; Q, Qe 2V )} (1.29) 
+ ( J d 6 zS [C 3 H/($, Q) + \f ab W aa yV*] + h.c.) . 

The invariance under Weyl symmetry is established by the following transformations 
for the compensator: 

C -)• e- 2r C. 

Because it does not appear in the action with derivatives, it is not a propagating 
degree of freedom and the action is completely equivalent to the case without the use 
of a compensator. This internal symmetry is broken when < C >^ (however the 
actual value does not matter as it can be chosen to be anything), but the breakdown 
is to nothing. 
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At the quantum level, however, the measure is not invariant due to the chiral 
anomaly [7j for matter fermions, 

<f[tt] -> d{^}exp{^ J dz 6 £ 2rW a W a + h.c, } (1.30) 

with 

c a = T(G a )-Y,T a (r), (1.31) 



r 



where T a (r) is the trace over the squared-generator matter representation a, and 
T(G a ) is the trace over the adjoint squared generator. In order to maintain local 
Weyl invariance, this anomaly must be canceled by the replacement 



3r 

/«(*) fa = /.(*) - ^ InC. (1.32) 



This ensures that the original SUGRA action of Equation (1.27) is equivalent to the 



Weyl symmetric action Equation (1.29). We can see that in the gauge C = 1, the 
effect is nil and we are in the original SUGRA frame. The Kahler -Einstein frame 
corresponds to the choice 

log C + log C = ^|harmonic , (1-33) 


where on the right we take only the piece that is the sum of chiral plus anti-chiral 
parts. 

We must also do matter field redefinitions in order for them to have canonical 
normalizations. To do this we expand the Kahler potential in matter fields 

if ($, $, Q, Qe 2V ) = K m ($, <5) + Z IJ -Q J e 2V Q I + .... (1.34) 
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For simplicity we consider a single matter field multiplet in representation r (the case 
of more fields is easily generalizable) and find that the kinetic terms are contained in 

J dz 8 ECCe- 1 3 Km Z r (^,^)Qe 2V Q, Z r = Z\. (1.35) 

If one does a field transformation of the form Q — > e TZ Q, the functional measure 
is again not invariant. This results in the Konishi anomaly [75] [HI] with another 
contribution to the measure of the form 



exp 



J dz 6 £ 2r z W a W a + h.c.y (1.36) 



This subsequently implies redefining the gauge coupling function: 



H. ^ /.(•) - InC - ^r z . (1.37) 



We then get canonical normalization of the matter kinetic term, 



J dz 8 ECCe-^ Km e TZ+fz Z r (^,^)Qe 2V Q, (1.38) 



by simply making the choice 



tz + tz = ln(CCe ^ Km Z r )\ harmonic . (1.39) 



This together with the appropriate gauge coupling redefinition, Equation (1.37), en- 
sures that the action remains locally Weyl invariant. 

The final step is to determine the gauge couplings and gaugino masses as they 
relate to the auxiliary components of all fields, viz., compensator, moduli, and mat- 
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ter fields. Denote the lowest component of the transformed gauge coupling function 



(Equation (1.37)) by h a = H a \, and its real part by h a R = $lH a \. The mass and gauge 



coupling are related by 



^ = M[F%f a (*)]\ - - ^F%(He^Z))\ (1.40) 



M[F%f a (<S>))\ - ^F%K m - ^F*ft(liiZ r ))|, 



where i is for all moduli of the underlying string theory. In going to the second line 
we have imposed the Einstein-Kahler gauge by using the F-component of Equation 



(1.33). The gauge coupling on the other hand is given by g a = h a R and in the 



Einstein-Kahler frame it is 



9a 2 = mm - ^~ 2 K m \ - ^lnZ r |. (1.41) 



Then, when considering F-type breaking, these expressions give the appropriate gaug- 
ino masses and gauge couplings after the Kahler potential and the original gauge 
coupling function (/ ) of the theory have been specified. 

Suppose the theory has a cut-off A, and we wish the renormalize the gaugino 
masses and the gauge couplings by evaluating them at a scale \i. This is done by 
shifting the Wilsonian gauge coupling function by a finite (1-loop) renormalization, 



H a ($, r, r z ) #„(<£>, r, Tz , y) = H a ($, r, r z ) - ^hi-, (1.42) 



and re-evaluating the expressions (here, the b a values are the (3- function coefficients). 
However, the gauge couplings we have derived until this point are not physical because 
the kinetic terms for the gauge fields are not normalized properly. We again must do 
a rotation into the proper frame, but this time of the gauge fields. The details are 
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omitted here, but the result is only to shift 3tH a by again using a new chiral superfield 
transformation parameters, Tz (see for details): 



$lH a -»• $H a - ^$tr v = (1.43) 



It is important to stress that this derivation of the soft mass contributions has 
made no impact on the scalar sector of the theory. The "Weyl" anomaly is often 
referred to as the "rescaling" anomaly referring to the fact that many authors have 
renormalized gauge couplings with an implicit compensator field associated with the 
mass scale of the theory, i.e., /3-function running appears usually as ln(^^) instead 



of as it is seen in Equation (1.42) (without \C\). It is this fact that has led authors to 
derive soft contributions for the scalars of the theory in addition to gaugino masses 
[94"] . However, it has been argued strongly in [46] that the problem with the usual 
derivation is that it is based on conformal invariance rather than strictly on Weyl 
invariance alone. It is ultimately argued in that paper that the Weyl anomaly does 
not contribute at all to scalar masses. 

In this thesis, a somewhat impartial approach is taken regarding the puzzles sur- 
rounding the anomaly. The models that are examined in the following chapters will 
have anomaly contributions as prescribed in the "usual" derivations (mAMSB and 
HCAMSB) that include scalar masses (resulting from anomaly rescaling) that evolve 
as 

< « -7?^3/2, (1-44) 



where 7, is the anomalous dimension. We also consider cases that contain anomaly 
soft contributions to gauginos as advocated by de Alwis [46] , as it is in the inoAMSB 



class. In all cases considered here, upon supersymmetry breaking, the relation (1.40) 
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implies soft gaugino masses of the for 



M a = ^m 3/2 , a = 1,2, 3, (1.45) 



with b a = (y , 1, —3). In any case, we will see that the various models will be distin- 
guishable at the LHC, at least at the lOOftT 1 level with y/s = 14 TeV. 



1.5 Computational Tools 

Renormalization Group Equations 

As described earlier in this Introduction, inputs at the high-scale (the GUT or string 
scales for instance) are required as boundary conditions for evolution of soft SUSY 
breaking parameters. These parameters then are used in determining physical mass 
parameters at the TeV scale by matching the MSSM renormalization group equations 
(RGEs) with low energy boundary conditions, i.e., weak-scale gauge and third gener- 
ation Yukawa couplings. 

All of the RGE parameter solutions in this research are obtained through an up- 
down iterative matching procedure using the Isasugra subprogram of Isajet v7.7S0 
|88j . Isasugra implements a full two-loop RGE running of the MSSM parameters 
in the DR-scheme using Runge-Kutta integration. The iteration starts with high- 
precision weak-scale gauge and Yukawa couplings, runs up to the the high scale where 
the user parameters are used as inputs, and then returns to the weak scale. At the 
end of the iteration, MSSM masses are recalculated and the RGE-improved 1-loop 
corrected Higgs potential is minimized. The /3-functions are also re-evaluated at each 

threshold crossing during each iteration. The iterations terminate at the prescribed 

8 This is derived explicitly for the inoAMSBcase in Se ction |4.2[ 

9 Isajet is publicly available code and can be found at http:77www.nhn.ou.edu/~isajet/ . 
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Model 


Parameters 


m AMQR- 
II1/\1V10I3. 


IilQ, ''t3/25 udJl Pj hoijliyjJ,) 


HCAMSB: 


a, 1713/2, tan (3, sign(fi) 


inoAMSB: 


m 3 / 2 ,tan (3, sign(fi), 




(M 1? M 2 ,M 3 ) = (6.6,1,-3) x m 3 / 2 



Table 1.3: 



level of convergence, which is that all RGE solutions except fi and Bfi are within 0.3% 
of the last iteration. Because of their rapid variation at the weak-scale, the latter are 
required to have convergence at the 5% level. 

Isajet v7.79 has the mAMSB and HCAMSB models coded into it. No extra coding 
was necessary for inoAMSB as the high-scale inputs are non-universal gaugino masses 
in the ratio 6.6:l:-3 as usual for AMSB, while all scalar and A-parameters inputs are 
highly suppressed, i.e., zero. For convenience, the parameters of each of the three 
models described in this work are listed in Table 11.31 

Event Simulation 

Isajet 7.79 - 7.80 is used for the simulation of signal and background events at the 
LHC A toy detector simulation is employed with calorimeter cell size At] x A<ft = 
0.05 x 0.05 and —5 < rj < 5. The hadronic calorimeter (HCAL) energy resolu- 
tion is taken to be 80%/\/E + 3% for \rj\ < 2.6 and forward calorimeter (FCAL) 
is 100%/v^ + 5% for \rf\ > 2.6. The electromagnetic (ECAL) energy resolution is 
assumed to be 3%/y/E + 0.5%. We use the UAl-like jet finding algorithm GET JET 
with jet cone size R = 0.4 and require that E T (jet) > 50 GeV and \rj(jet)\ < 3.0. 
Leptons are considered isolated if they have pr(e or fi) > 20 GeV and \r]\ < 2.5 with 
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visible activity within a cone of AR < 0.2 of TjE^? Us < 5 GeV. The strict isolation 
criterion helps reduce multi-lepton backgrounds from heavy quark (cc and bb) pro- 
duction. 

We identify a hadronic cluster with Et > 50 GeV and \rj(j)\ < 1.5 as a 6-jet if it 
contains a B hadron with Pt(B) > 15 GeV and |?7(-B)| < 3 within a cone of AR < 0.5 
about the jet axis. We adopt a 6-jet tagging efficiency of 60%, and assume that light 
quark and gluon jets can be mis-tagged as 6-jets with a probability 1/150 for E T < 100 
GeV, 1/50 for E T > 250 GeV, with a linear interpolation for 100 GeV< E T < 250 
GeV [65]. 

Isajet is capable of generating events for a wide variety of models. Once the pa- 
rameters of the theory are defined, RGE evolution determines weak scale masses and 
2 — > 2 processes are weighted by their cross sections and generated. 

In addition, background events are generated using Isajet for QCD jet production 
(jet-types include g, u, d, s, c and b quarks) over five pr ranges as, for example, 



in Tables 3.4 and 4.3 Additional jets are generated via parton showering from the 
initial and final state hard scattering subprocesses. Also generated are backgrounds 
in the W + jets, Z + jets, tt(173.1) and WW, WZ, ZZ channels. The W+ jets and 
Z + jets backgrounds use exact matrix elements for one parton emission, but rely on 
the parton shower for subsequent emissions. 

Dark Matter 

For all the models considered in this work, to calculate the relic density of neutralinos 
and direct detection rates we used, respectively, the IsaRED and IsaRES subroutines 
of the IsaTools package found in Isajet (55] . To calculate indirect detection rates 
following from neutralino scattering, annihilation, and co-annihilations, the Dark- 
SUSY 10 [60J package was used. DarkSUSY, by default, uses the Isasugra subprogram 



10 Publicly available at http://www.physto.se/~edsjo/darksusy/ 
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of Isajet v7.69 [EE] to calculate the SUSY spectrum, with exception to Higgs masses 
for which it uses FeynHiggs [61J. In order to perform the calculations for our particular 
models, we transplanted the Isajet v7.79 code into the DarkSUSY package. 
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Minimal AMSB 



With the origin of the anomaly contributions to the particle spectrum understood we 
can now look at models with high scale inputs that are connected to the TeV scale 
through renormalization group running. The first model we encounter has been ana- 
lyzed extensively in the literature [30] [87] and will serve as the model of comparison 
for the other AMSB models to be considered. The model is known as the Minimal 
Anomaly Mediated Supersymmetry Breaking Model (mAMSB) and its most impor- 
tant aspects, relevant for LHC and DM considerations, are highlighted in this chapter. 

The mAMSBmodel has the attractive feature that it depends on only a few GUT- 
scale input parameters: 

m 0i m 3/2, tan (3, and sign(/i), (2.1) 
each of which have been mentioned in Chapter [TJexcept for m , which will be discussed 



shortly. As discussed in Section 1.4, the AMSB contributions to the gaugino masses 
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are given by 



33 Jh_ 2 (2 2) 

5 16vr 2 3/2 V ; 

M 2 = A- 2 m\ l2 (2.3) 

M 3 = -3^*4/2, (2-4) 



where <?j are the running MSSM gauge couplings. Furthermore, the scalar masses and 
trilinear parameters are given by 



m ) = ~\{jgPa + ^^/} m 3/2 and ( 2 - 5 ) 

/ 167T 2 



where /3 g and /3/ are, respectively, the gauge coupling and Yukawa coupling /3-functions, 
and their correspond anomalous dimensions are denoted by 7. Note that the above 
soft terms are parametrically tied to supersymmetry-breaking through m^/2 (Section 



1.4). Also note that all of the above equations are valid at any scale. 

There are several phenomenologically important points to be made about the 
AMSB contributions. 

i. Anomaly contributions to scalars are the same for particles with the same quan- 
tum numbers, while first and second generation Yukawa couplings are negligi- 
ble. In the case that AMSB dominates the soft contributions, flavor-violation 



is safely avoided. In Figure 2.1 the running of the soft parameters Mi, M2, M3 



the third generation squark doublet mass {^ hrfk ), the right-handed sbottom 



mass {y m f ), and the bottom trilinear parameter (Ab) are shown. There is 
implicit dependence on Ab in both of the squark mass parameters. 

ii. The anomaly contributions to scalars are partially determined by anomalous 
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dimensions, which are negative for sleptons. After RGE running to the weak 



scale, sleptons remain with negative masses (see Figure L3 in 2.1 ). The minimal 



solution to this problem is to introduce the ad hoc parameter, mo, in 2A_ at the 
GLT-scale to prevent tachyonic sleptons at the weak scale. This is achieved 
simply by adding m\ to all scalar squared-masses and the mo value is taken as 
a free parameter of the mAMSBmodel. 

Hi. Scalars with different helicities are nearly (left/right) degenerate at the weak 
scale in the mAMSBmodel. This is despite possible splitting at the GUT scale. 



This is demonstrated in Figure |2.2| where, for example, selectrons and up- 
squarks masses are run from the GUT scale to the weak scale for a mAMSB 
point with m =300 GeV, m 3 / 2 = 50 TeV, tan = 10, and fi > 0. While left 
and right sparticle masses are split at the high scale, the evolution (accidentally) 
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drives them closer at lower scales. 
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Figure 2.2: Demonstration of near-degeneracy of left and right superparticles at weak 
scale in the mAMSB scenario for the point (m , m 3 / 2 , tan 0) = (300 GeV, 50 TeV, 
10) with fj, > 0. 



iv. Requiring weak-scale gaugino masses places m 3 / 2 at much higher values than 
in other supersymmetry-breaking scenarios. With 777.3/2 ~ 0(10 — 100) TeV 
range, gravitinos have a short enough lifetime to avoid the gravitino problem 



v. After evaluating gaugino masses at the weak scale (Equations (2.2 - 2.4) 



|M 2 | is invariably the lightest of them (see Figure 2.1). This, along with the 



potential minimization conditions, affects the diagonalization of the neutralino 



and chargino mass matrices of Equations (1.19) and (1.22) to produce uwo-like 



lightest mass-eigenstates in both cases. Thus Z\ and W% have nearly degenerate 
masses in AMSB type models. 



34 



This last point is particularly important because it potentially allows for a dis- 
crimination between AMSB and other types of supersymmetry-breaking effects at 
the LHC. Denote the mass difference between the lightest chargino and the lightest 
neutralino, — , by AM| . As AMg decreases, the chargino decay width 
decreases and its lifetime becomes longer. Thus, there is a possibility of detecting 
W±s as highly ionizing tracks (HITS) in LHC detectors. 

AM^ ~ O(MeV) splittings are rather small but W\ can usually decay to e or 
(jl through the three-body process W\ —> evZ\. When the gap opens beyond the 
pion mass (~ 140 MeV) hadronic decays of the chargino are allowed [32] [78]. Further 
enlarging the mass gap leads to channels with 2n and 3tc final states. The W\ width 
would be smallest when it decays to a single lepton and obviously grows as other 
channels open. Chargino detection then falls into two categories [30J: 

• AM% < 771^+ : 7T decay modes are unavailable and W\ can leave a track far out 
in the muon chambers! 

• m n + < AM % < 200 MeV : high charginos decay in the inner detector region. 
The decays W x -> vr ± Z 1 and W x -> lvZ x are accessible with I and tt emitted 
softly. This regime has a clear signal. The charginos appear as highly ionizing 
tracks in the inner detector and decays to SM particles. But the SM particles 
are too soft for energy deposits in the Ecal or the Heal and appear to contribute 
to the missing energy already taken away by the Z\. This leads to an observable 
HIT that terminates, or a stub, and large amounts of p T . 

We will come back to this in Chapters [3] and [4] when we discuss the LHC phenomenol- 
ogy of AMSB models at greater depth. 
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Hypercharged AMSB 



3.1 Introduction to the HCAMSB Model 

Hypercharged anomaly-mediation is composed of two mechanisms that induce masses 
for visible sector matter fields. The soft masses come from the anomaly mediation al- 
ready discussed and an additional U(l) mediation, and the latter depends intricately 
on a particular D-brane setup. While this form of mediation is not general in the 
way that anomaly contributions are, it is an interesting pathway to understanding 
how D-brane constructions can impact the visible sector. In the next section the 
description of how the U (1) gives masses to MSSM fields is given. Then the U(l) will 
be paired with AMSB to give the full HCAMSB model. Following the next section 
the phenomenology of the model will be discussed. 
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3.2 Geometrical Setup with D-branes 

There should be some comments made about D-branes first because the model re- 
lies on them. D-branes are extended objects on which strings can terminate with 
'D'irichlet boundary conditions. For this discussion we consider only type IIB string 
theory with Dp-brane^j] that have p = 1 — 9, odd. These branes fill the usual 3+1 
space-time, but can have superfluous dimensions that must be wrapped by internal 
cycles to make them effectively invisible. A D7 brane, for example, requires a 4-cycle 
wrapped within the internal geometry, and D5 requires a 2-cycle, etc. 

D-branes have interesting features that are useful for constructing realistic field 
theories. Most importantly for this discussion are the following two properties: 1.) 
chiral matter can exist as open strings at the intersection of two D-branes, and 2.) 
interactions with local curved geometry lead to bound states of D-branes known as 
"fractional branes" when branes of different dimensions are involved. In this model 
the bound states of branes occur at singularities in the Calabi-Yau (CY) manifold. 
There are two CY singularities; the brane located at one of the singularities will be 
the visible brane, and the other will be hidden. In addition, the two branes will share 
properties that allow for the U(l) mediation. 

U(l) Mediation of SUSY Breaking 

The main idea behind U(l) mediation is that, given a proper geometrical setup, a 

brane can communicate SUSY breaking to another brane despite there not being any 

open string modes connecting them. An F-type breaking occurring on one brane can 

be mediated to another through bulk closed-string modes that have special couplings 

to gauge fields on the branes. To understand the mechanism, first consider a single 
1 Dp-branes have d-p-1 Dirichlet and p+1 Neumann boundary conditions. 
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D5 brane separated within the CY manifold with a U(l) symmetry and associated 
gauge field A. Dp-branes are themselves sources of generalized gauge fields C p+ i, 
for p=l, 3, & 5 in IIB theories [50] and exist in the bulk. C p+ i have induced linear 
Chern-Simons couplings to the brane gauge fields. In particular, C4 couples to A 
through 

C cs = C 4 AdA + AAF 5} (3.1) 

where F 5 is the five-form field strength of C4. F 5 has the special property of self-duality 
in 10 dimensions, that is, F 5 = *F 5 . The equations of motion for this Lagrangian are 

dF 5 = dA A S brane , and * F 5 = dC 4 + A A S brane , (3.2) 

Now consider a 2-cycle a wrapped by the D5 brane, and the four-cycle (3 dual to 
a. When the extra dimensions are reduced, C4 leads to a massless two-form for each 
two-cycle wrapped by the brane, and each 4-cycle leads to a massless scalar. For the 
cycles a and (5 this means that we have 

C = j C 4 (2-form) (3.3) 

J a 

<p = / C4 (scalar), (3.4) 



and we assume that there are no other cycles around. C and ip are related via the 
self-duality of F 5 . A unique basis can be chosen for the expansion of F 5 consisting of 
a 2-form a and its dual 4-form p that satisfy the following properties: 



a = l, / p — 1, and / a A p = 1. (3.5) 

a J p JCY 
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These forms are related by Hodge duality: 



* 6 p = p 2 a, [ a A * 6 a = — , (3.6) 

JCY ^ 



where p, ~ string scale and characterizes the size of the compactification. We can 
expand F5 in this basis and KK reduce it: 

F 5 = dC A a + (d(j) + A) A p + .... 

From this self-duality of F 5 is satisfied provided that 

H dC = p 2 {d<f) + A), (3.7) 

where p 2 appears in the a and p duality relations. This equation is a solution to the 



equations of motion of the action dual to the 4D version of Equation (3.1 ) + C-kinetic 
term. This results in a low energy mass term for A. 

Now when we consider the actual setup which is similar to that already considered 
but consists of two D5 branes, one visible (V) and one hidden (H). There are U(l) 
gauge bosons on the branes denoted now by Ay and Ah, and each D5 brane wraps 
its own two- and four-cycles (and associated 2- and 4- forms), denoted by ay, an, 
j5y, and (3h- Actually we choose the CY geometry such that these cycles are topolog- 
ically the same. For instance, if ay and an are topologically the same two-cycle they 
can be continuously transformed into each other. If we follow the same procedure 



outlined above, we again arrive at Equation |3.7| but with A — > Ay + Ah, and this 
leads to a mass for the combination Ay + Ah- With string scale compactifications, 
this combination is quite heavy, and is lifted from the low-energy spectrum. However, 
the remaining light combination, Ay — Ah, does survive to low energies as a light 
vector boson [IB]. In this model, this combination is identified with the ordinary 



hypercharge boson, and the effects of supersymmetry breaking are imparted to the 
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superpartner of Ah- Thus, the bino acquires extra soft contributions that will alter 
the usual mAMSBcontributions in interesting ways to be described in the next section. 



3.3 Spectrum and Parameter Space 

In this section we examine the mass spectrum and understand how it evolves from 
GUT-scale running to the weak scale. Because we have the Feynman rules for the 
MSSM, all rates can be calculated as long as the masses are given. This includes 
rates of rare processes that can be sensitive to the choice of model parameters. These 
measurements are used to place constraints on the parameter space of the theory. We 
will also constrain the parameter space as much as possible from other theoretical 
considerations such as the requirement of proper electroweak symmetry breaking. 
The soft mass contribution RGEs are 



33 n z 

Mi = M 1 + Y ^- 2 mi /2 (3i 



M a =^mi /2 a = 2,3 (3.9) 

^ = \^ + Tf^ m ^ (3 ' 10) 

A f = ^m 3/2 . (3.11) 



The difference between these and mAMSB renormalizations is that the equation for 
Mi has an extra Mi input at the high scale (GUT scale chosen for convenience) that 
accounts for the hypercharge contribution. Winos and gluinos only receive negligible 
two-loop contributions from the bino [38]. To compare the U(l) effects with those of 
AMSB we rewrite Mi proportional to m 3 / 2 as 

Mi = a m 3/ 2 
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so that the bino soft mass reads as 

33 a 2 

M 1 = (a + -^- 2 )m 3/2 . (3.12) 

The parameter space of the HCAMSB model is then 

a,m 3 /2,tan (3, and sign{[i), 

which resembles the mAMSB space except that a will be the parameter that helps 
to avoid tachyonic sleptons at the weak scale. 

Figure [3~1~1 shows the particle mass spectrum at the weak scale for both the mAMS- 
Band HCAMSBmodels against m and a respectively. In order to compare the effect 
of mo with that of a (= we consider its ratio to m 3 / 2 , while fixing 7713/2 = 50 

TeV and tan /3 = 10. We will fully explore the impact of varying the other parameters 
later. 

These plots have a few qualitative similarities. The first feature is that each 
has a yellow-shaded region corresponding to RGE solutions with tachyonic sleptons. 
These regions are forbidden because the scalar potential should not be minimized 
by charged scalars, as this would lead to breakdown of electric charge conservation. 
These regions are where they would be expected; for mAMSBthis is near m ~ 0, 
and for HCAMSBit is around a ~ where the bino contribution is small and pure 
AMSB is recovered. In general, as each parameter increases the general trend is for 
masses to increase (although here are important exceptions in each case). Each has 
nearly-degenerate Z\ and W\ and their masses remain relatively flat for all mo and 
a. This is also the case for the g. Conversely /x is seen to decrease with increasing 
values of the parameters. The upper edges themselves are due to improper breaking 
of the electroweak symmetry which is signaled by /i 2 < 0. 

There are also notable differences between the spectra that eventually lead to dis- 
tinct phenomenology. In the mAMSBcase, there is a near-degeneracy between left 
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-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 

a 

(b) 

Figure 3.1: Mass spectrum for mAMSB(left) and HCAMSB(right) with model pa- 
rameters m 3 / 2 = 50 TeV, tan (3 = 10, and \i > 0. 
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and right particles of the same flavor due to the nearly-equal /3-functions and that 
is universally added to all scalar squared masses. HCAMSBon the contrary exhibits 
a left-right split spectrum. For example, mAMSBhas mg R — m^ L but these values 
can be seen to differ by over 0.3 TeV in the HCAMSBcase for all a. 

The plot also shows that stop and sbottom masses actually decrease in HCAMS- 
Bfor larger hypercharge contributions unlike the other scalars of the theory. Within 
mAMSBall scalars increase with ml because this contribution is simply added to all 
high-scale, scalar, squared soft mass values . 

The parameters /i, M\, and M 2 determine the composition of the neutralinos 
and are different between mAMSBand HCAMSBmodels. These parameters mix to 
form the eigenstates of the neutralino mass matrix, Zi. Because the values of these 
parameters and their relative ordering determine the composition of the LSP, they 



are also responsible for its interaction properties. For example, Table 3.1 shows the 
ordering of neutralino mass parameters and the main components of the neutralino 
mass eigenstates for the benchmark point for mAMSBand Point 1 for HCAMSB(the 
selection of representative points will discussed in the next section). Since M 2 has 
the smallest value in each case the LSP is wino-like with a small mixture of bino and 
higgsino components. However, Z 2 is bino-like for mAMSBand mainly higgsino for 
HCAMSB. We can see already that when Z 2 is produced in collisions that its decays 
will be heavily model-dependent and will lead to final states with either strong bino 
or higgsino couplings. This will be crucial is distinguishing between the two models 
in the LHC section. 

We should now turn our attention to understanding the sources of the spectrum 
patterns described in the previous few paragraphs. In order to understand mass pa- 
rameters at the weak-scale we need to examine how they evolve from the GUT scale. 
It is seen from the RGEs that the AMSB contribution to scalar masses are determined 
by anomalous dimensions, 7^. Because ji are positive for squarks and negative for 
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Neutralino 
Eigenstates 


mAMSB 


HCAMSB 


fi> Mi > M 2 


M l > fi > M 2 


z 2 

Z* 

z A 


99.09 % wino; ~ 1 % bino-higgsino 
99.35 7o bmo; < .7 7o higgsmo-wmo 

99.72 % higgsino; < .30 % bino-wino 

98.73 % higgsino; < .30 % wino-bino 


99.04 % wino; ~ 1 % bino-higgsino 

70.30 % higgsino; 29.09 % bmo 
99.74 % higgsino; < .30 % bino-wino 
70.88 % bino; 28.93 % higgsino 



Table 3.1: Compositions of the neutralino mass eigenstates. Both models have wino- 
like LSP, but heavier states differ due to order of //, M%, and M 2 . 



sleptons, the masses of particles begin at the GUT scale with their respective signs. 
Each scalar receives a contribution from the hypercharge mediation of the form jl8] 



*»?«) - -^ m HjL) (3 - i3) 



where is the hypercharge of the i th scalar. This contribution serves to uplift masses 
as the scale Q decreases. This accounts for the large left-light splitting between the 
scalars in figure 3.1[ since right-handed particles have a greater hypercharge in gen- 



eral. As Q approaches the weak-scale the masses become larger and Yukawa terms 



tend to dominate the RGEs which leads to suppression of masses. Figure 3^ shows 
these effects for hypercharge anomaly mediation and for pure AMSB. In the case of 
pure AMSB, the sleptons begin with negative mass parameters and are tachyonic at 
the weak scale as expected. Sleptons have large hypercharge absolute values, and the 
figure shows that the hypercharge contribution lifts slepton masses to positive values 
for the HCAMSB case. 

Pure hypercharge mediation however has the opposite effect on the third gener- 
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Figure 3.2: HCAMSB (solid) and pure AMSB (dotted) mass evolution from GUT 
scale to the weak scale for third-generation sleptons, third-generation squarks, and 
the up-type Higgs. 

ation squarks. As we evolve to the weak scale, the left-handed stops may become 
tachyonic because they have relatively small hypercharge values and the same large 
Yukawa couplings of top quarks. But the stops also receive AMSB contributions that 
leave them with positive masses in the HCAMSB scenario. The anomaly and hy- 
percharge mediation mechanisms are complementary in the respect that each helps 
to avoid the tachyonic particles that are present when only the other mechanism is 
present. 

We can understand the extra suppression with increasing a for third generation 
squarks by taking a closer look at the RGEs. Both ti and bi belong to the same 
577(2) doublet, Q%. The running of the doublet mass includes the terms [20] 




3 Y t 2 X t + Y b 2 X b 



(3.14) 



d{log Q) 
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where X tj b is defined to be 

X t = rr4 + m\ R + m 2 Hu + A 2 (3.15) 
X b = m 2 Q3 +ml R + m 2 Hd + A 2 t . (3.16) 

Again, right-handed particles have larger hypercharge than SU(2) doublets, and the 
values m 2 R and m~ steepen the running slope near the weak scale. Then it is evident 
that for larger values of a in HCAMSB, the third generation doublet receives extra 
suppression from Yukawa effects relative to other generations. This suppression affects 
tir and ^-production rates at the LHC for the HCAMSBmodel as will be discussed 
in the next section. 

We can also see which RGE effect leads to /i-suppression with increasing a. The 



tree-level scalar minimization given by Equation 1.17 shows that /i goes as —m 2 H u 
and increasing m 2 H implies decreasing /j. The m 2 H RGE includes the terms 



ii^ 5 -h- M ' +3Y ' x '- < 3 - 17 > 



Again the large Mi value uplifts m 2 H in the early running from Q = Mqut, and at 
low Q Yukawa effects again dominate over the hypercharge effects. This is shown in 
Figures |3.2 and 3J3 In Figure 3^3 three values of a have been chosen and the 1-loop 



RGEs have been used in the evolution. As a increases the Higgs mass increases at 
the weak scale. For the largest a shown the Higgs mass is positive which could imply 
negative fi 2 and therefore no electroweak symmetry breaking. However, when large 
1-loop corrections are added to the scalar potential /1 2 is once again positive. In Fig- 



ure 



3.3[ a ~ 0.195 is the upper edge of parameter space beyond which EW-breaking 
does not occur. 

As already mentioned, /j is an important parameter in the neutralino sector. To- 
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wards higher values of a, fi decreases as mentioned above and moves nearer to the 
value of M 2 (over all of the parameter space Z\ is mainly wino and M% ~ M 2 ). 
Here the Z\ mass state is a mixture of wino and higgsino. Because of this, the mass 
splitting AM% between W\ and Z\ will increase leading to a shorter lifetime of the 
former. 

To understand this last point, note that in the limit of /i, M\ » M 2 the mass 
eigenstates Z% and form an SU(2) triplet with common mass M 2 . The symme- 
try is broken by gaugino-higgsino mixing which leads to mass splitting between the 
neutralino and the chargino. As we saw in Chapter |2j AM^ is important in the 
detection of charginos. We now see that the detection of the chargino can depend on 



the value of a, and its lifetime is shown in Figure 3^ As can be seen, larger a leads 
to shorter W\ lifetimes. 



All of these results can be summarized as in Table 3.2 for the mAMSBpoint with 
(mo, m 3 / 2 , tan (3) = (300 GeV, 50 TeV, 10), and for similar parameter choices for two 
HCAMSBpoints, with a low and high value of a, 0.025 and 0.195, defined respectively 
as points HCAMSB1 and HCAMSB2. By comparing the model lines we see that the 
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Figure 3.4: W\ lifetime as a function of a. 

table reflects all of the HCAMSBfeatures: L-R splitting, wino-like Z-y/Wi, mq 3 - and 
/^-suppression and AM^ increase with increasing a, etc. 



3.4 HCAMSB Model Constraints 

We now examine the allowed parameter space of the HCAMSBmodel and begin to in- 
fer some sparticle mass ranges observable at the LHC There are both theoretical and 
experimental constraints for the model parameter space. On the theory side we have 
the requirement of proper electroweak symmetry breaking which implies that only 
m#2 can become less than zero because the scalar potential cannot be minimized by 
fields with conserved charges. If other mass parameters (including //) become negative 
and/or H u does not, then the corresponding region of parameter space is prohibited. 
Experimental constraints come in two varieties: direct and indirect. The direct con- 
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parameter 


mAMSB 


HCAMSB1 


HCAMSB2 


a 


— 


0.025 


0.195 


m 


300 


— 


— 


m 3 / 2 


50 TeV 


50 TeV 


50 TeV 


tan (5 


10 


10 


10 


Mi 


460.3 


997.7 


4710.5 


M 2 


140.0 


139.5 


137.5 


/i 


872.8 


841.8 


178.8 


m§ 

y 


1109.2 


1107.6 


1154.2 




1078.2 


1041.3 


1199.1 




1086.2 


1160.3 


2826.3 




774.9 


840.9 


427.7 




985.3 


983.3 


2332.5 




944.4 


902.6 


409.0 




1076.7 


1065.7 


1650.7 




226.9 


326.3 


1973.1 




204.6 


732.3 


3964.9 


vv 2 


879.2 


849.4 


233.1 




143.9 


143.5 


107.1 


m z 4 


878.7 


993.7 


4727.2 


m z 3 


875.3 


845.5 


228.7 


m z 2 


451.1 


839.2 


188.6 


m z, 


143.7 


143.3 


105.0 


m A 


878.1 


879.6 


1875.1 


m h 


113.8 


113.4 


112.1 



Table 3.2: Parameters and masses in GeV units for three case study points mAMSB, 
HCAMSB1 and HCAMSB2 using Isajet 7.79 with m t = 172.6 GeV and n > 0. 



straints considered are the LEP2 mass limits on or m^. Indirect constraints come 
from measurements that would be sensitive to supersymmetric particles appearing in 
loops. The indirect constraints considered here are the branching ratio from inclusive 
radiative B-meson decays, BF(b — > 37) [M], and (g — 2) M measurements. Consider- 
ation of cosmological constraints are postponed until Chapter [5] where Dark Matter 
in AMSB models is discussed. 

We begin exploring the parameter space by plotting various masses in the m 3 / 2 — a 
plane in Figure 3.5| The plots in the figure have excluded regions for improper EW 
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tan(P) = 10 tan(P) = 40 




-0.2 -0.1 0.1 0.2 -0.2 -0.1 0.1 0.2 

a a 



Figure 3.5: Allowed m3/2-« parameter space for tan (3 — 10 and 40 with /i > and 
m t = 172.6 GeV. The orange vertical lines at lower values of 7713/2 represent the LEP2 
excluded region where < 91.9 GeV. The yellow regions at the far-right, far- left, 
and center are regions where electroweak symmetry is improperly broken. The white 
regions are acceptable, and constant-mass contours for Ui,g,ti, and are shown. 



symmetry breaking (yellow-thatched region) due to tachyonic sleptons at a ~ and 

fi 2 < at the extreme a values. The orange region at lower m 3 / 2 values is where 

is below the LEP2 limit of 91.9 GeV [80] in the search for nearly degenerate V^is and 

The plot shows that the region with 777.3/2 < 30 TeV are excluded by the LEP2 

chargino limit. For m 3 / 2 ~ 30 TeV, we have for the gluino m g ~ 730 GeV which is 

2 The LEP2 limit on the Higgs mass is mH {SM) > UiAGeV. While this limit is possibly con- 
straining there is an estimated ±3 error in the calculation of m^. Since > lllGel^ over the 
entire allowed parameter space we do not show this constraint in Figure 3.5 
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a) m 3/2 = 50 TeV, tanp = 10 
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c) m 3/2 = 100 TeV, tanp = 10 
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b) m 3/2 = 50 TeV, tanp = 40 
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d) m 3/2 = 100 TeV, tanP = 40 
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Figure 3.6: The branching fraction for b — > 57 vs. a for combinations of 7713/2 = 50 
and 100 TeV and tan /? = 10 and 40. 



beyond the reach of any reasonable search at the Tevatron, so the discovery potential 
for this model must be investigated for the LHC There is not an upper-bound on 
the gravitino mass, but the plot extends up to m 3 / 2 ~ 150 and 160 TeV for the 3 
TeV contours in rrig and rriu L respectively. These are somewhat above the reach of 
2.1 TeV for squarks and 2.8 TeV for gluinos [30] predicted in the case of mAMSBfor 



the LHC. We will explore the reach of this model in Section 3.6 



We also check whether there are regions of parameter space that agree with indi- 



rect measurements. Figure 3^3 shows the branching fraction of b — > 57 as a function 
of a for four pairings of the m.3/2 and tan /3 parameters to span the space. The dashed 
red line represent the HCAMSBvalues calculated with using the Isatools subroutine 
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ISABSG [IB]. The blue line is the theory result for the SM at order a 2 s [M] with the 
range of values BF(b — > 57) = 3.15 ± 0.23 x 10~ 4 . The black-dotted lines are the 
combined experimental branching fraction results of CLEO, BELLE, and BABAR 
[29] and have values in the range BF(b — > S7) = 3.55 ± 0.26 x 10 -4 . The plot shows 
that there are regions of near-agreement in the parameter space of each case. Excep- 
tions include the low a region in frame b) where the BF is too big, and very high a 
values in frames a), c), and d) where BF is too small. Finally, we close this section 
with contribution to < - 3 ~ 2 ^ M , denoted as Aa^ USY , in the HCAMSBmodel. Figure 



3.7 



shows Aa^ USY calculated with ISAAMU from Isatools [IT], with four plots with high 
and low values of 7713/2 and tan /3. In each of the four cases, low a leads to relatively 
light fiL and vl that appear in loop corrections to the photon vertex. This appears 
in the plots as larger corrections at low a to SM predictions. Parameter values with 
Acl susy > go x 10~ 10 or < are disfavored [20]. 



3.5 HCAMSB Cascade Decays Patterns 

All of the signatures at the LHC for the HCAMSB model will be shown in this sec- 
tion. That is, the physical outcome of those regions of parameter space that are not 



already excluded will be explored. We will also rely on the findings of section 3J3 to 
understand the production rates when necessary. 

We begin by examining the largest LHC production cross sections in order to 
understand what are the HCAMSBsignatures. Note that all of the following analysis 



is done for the 14 TeV pp beams. Table 3.3 shows for the three model lines the 
total SUSY production cross section and the percentages for the pair-production for 
gluinos, squarks, EW-inos, sleptons, and light stops. It is seen that the EW-ino pairs 
dominate over all other production rates while squarks and gluinos are produced in 
lower, but still significant, amounts. 
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m 3/2 = 50 TeV, tan(P) = 10 
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Figure 3.7: The SUSY contribution to (g — 2) M as a function of a. 

The dominant EW-ino cross sections come from the reactions pp — > W^W^X and 
pp — > WfZ\X. However, the W\ decays and the Z\ do not lead to calorimeter signals 
that can serve as triggers for LHC detectors. So we must instead look at squark and 
gluino production mechanisms. 

We first consider that 7713/2 = 30 TeV, the lowest value not excluded by exper- 
iment, and qualitatively discuss the emergence of the final states. Because of the 
absence of third-generation partons in the initial state, squark and gluino rates are 
determined by SUSY QCD and only depend on their respective masse^J Since m^, 
rriq R and rrig have similar values for low a, the final states gg, gq, and qq are produced 

at similar rates. In general, R-squarks are heavier than their L-partners due to the 
3 For production/decay rates see Appendices of Ref. [20] , 
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mAMSB 


HCAMSB1 


HCAMSB2 


a[fb] 


7.7 x 10 3 


7.4 x 10 3 


1.8 x 10 4 


g, q pairs 


15.0% 


15.5% 


14.3% 


EW-ino pairs 


79.7% 


81.9% 


85% 


slep. pairs 


3.7% 


0.8% 




Mi 


n AO/ 
(J. 47o 


0.Z7O 


r [-0/ 

5.5% 


BF(Z 2 Z X Z) 


0.01% 


7.7% 


22.3% 



Table 3.3: HCAMSBrates for the LHC. 



U(l) contribution. Already for low a, ra^ R and m £fl > rrig (similarly to m 3 / 2 



50 



TeV in Figure 3.1). The subsequent decays of these squarks enhance the production 
of gluinos through ur — > ug and cr — > eg. Gluinos finally decay in quark-squark 
pairs, and they have the highest rates into bb\ + h.c. and tti + h.c. and subdominant 
rates into other qq^ pairs. 

Conversely, as a increases, right-handed sparticle masses become larger to the 
point that eventually they cannot be produced in collisions. At higher values of a, 
left-sparticles become heavier than g, while t± and h\ are significantly lighter and 
again can be found in the main quark-squark decay modes of gluinos. At the highest 
a values, gluinos decay only into quark-squark pairs involving tis or &is: g — > tt\ or 

We also find for high a that, in addition to g — > ti/&i + X, direct production of 
b\bi + h.c. and trfi + h.c. pairs dominate over gg,gui, and gdi- The ti and b\ mass 
eigenstates are mainly L-squarks at high a, and appear approximately in a weak dou- 
blet. Thus they have nearly the same mass, ~ mg i , and their production rates 



are nearly identical. Figure 3.8 shows the direct production cross section of sbottom 
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LHC: pp collisions, V7= 14 TeV 
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Figure 3.8: Total b\b\ production cross section. 

pairs using the program Prospino2.1 [33], and it is understood that pair-production 
cross section of stops is nearly equal. 

To recap the findings of the previous three paragraphs we see that light stop and 
sbottom squark, as well as top and bottom quarks are produced in the following ways: 

low a: • gluino production and main decays to quark- squark pairs; 
high a: • gluino production and subsequent decay purely to quark-squark pairs; 
• direct stop and sbottom pair production. 

Obviously the production of stops and sbottoms is important in HCAMSBphenomenol- 
ogy. Then to proceed we need to examine the decay patterns of these particles to 
arrive at the final state. 

We now move from 7713/2 = 30 TeV up to 50 TeV to match the benchmark points 
and we find that the features of the previous paragraphs are unchanged. To proceed, 
the branching fractions of the lightest sbottoms and stops are plotted as a function 
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Figure 3.9: Branching fraction of h\ for model parameters m 3 / 2 = 50 TeV, tan 
and sgn(fi) > 0. 
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Figure 3.10: Branching fraction of ti for model parameters m 3 / 2 = 50 TeV, tan 
and sgn(fi) > 0. 
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of a in Figures 3J3 and 3.10 The features of low a are simple: both light sbottoms 
and stops decay to lightest charginos with about 67% and to lightest neutralinos 
with about 33%, and those make up the entire branching fractions. As \a\ increases, 
ft decreases and approaches M 2 , and again Z\ makes the transition from wino-like 



to a wino-higgsino mix. Z 2 and Z 3 are mainly higgsino (as in Table 3.1) and be- 
come lighter with decreasing /1, so that the decay channels b± — > bZ 2 ,bZ 3 ,tW 2 and 
ti — )■ tZ 2 ,tZ 3 ,bW 2 open. Finally, at the largest values of \a\, the sbottoms and stops 
decrease in mass and the decay modes close. In particular, it can be seen in Figure 



3.10 that all two-body modes for the light stop are closed, but its decays proceed via 



the t\ — > blvZ\ and cZ\. 

We also check what effects increasing m 3 / 2 has on the branching fractions. Figure 



3.11 shows the branching fractions verses m 3 / 2 for sbottoms and stops at low |a|. It is 
seen in both cases that out to very large values of m 3 / 2) two-body decays modes with 
Wi and Z\ are dominant while all others are subdominant. Explicitly, these modes 
are t x -» bW i ,tZ 1 and 61 ->■ bZ 1 ,tW 1 . 

Then, to summarize the HCAMSBsignatures, we have seen that gluino and squark 
production always lead to cascade decays to third-generation quarks and squarks. 
From these products should emerge high multiplicities of b — jets as well isolated 
leptons emerging from the t — > bW decays. There should also be significant $ T due 
to the presence of Z\S and us escaping the detectors. In addition to these, there is 
also the possibility of detecting highly-ionizing tracks as is usually for AMSB mod- 
els. This is due to the wino-like nature (M 2 << Mi, ft) of Wi and Z\ states and to 
their near mass degeneracy that leads to a narrower width for the former. Because 
Wi — > n + Zi has too soft a pion to leave a calorimeter signal (Chapter [2k, the HITs 
should be seen to abruptly terminate without a calorimeter signal. 
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Figure 3.11: Branching fraction for h\ (a) and t\ (b) vs. m 3 / 2 for model parameters 
a = 0.025, tan (3 = 10, and sgn^fi) > 0. 
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3.6 HCAMSB and the LHC 

In this section we analyze events using the Isajet 7.79 [BE] event generator (for specific 



details of the simulation, see section 1.5 of the Introduction). As we saw in section 



3.4[ the lowest gluino mass is too large to be detected at the Tevatron; we then must 
understand how HCAMSBwould appear at the LHC. It is expected that the LHC 
phenomenology will have similarities to mAMSB on general grounds. So part of the 
aim of this section will be to see in which ways this theory will be distinguishable 
from mAMSB at the LHC. 

The discussion in this section is split into two subsections. In the first we flesh 
out the final states of cascade decays and optimize the cuts used to observe them. In 
the subsection that follows we analyze the reach for the LHC. 

Final State Analysis 



We generate 2M events for the points HCAMSB 1 and HCAMSB2 from Table [372 
and compare with SM backgrounds. QCD jets are generated with Isajet and include 
the following jet types: g, u, d, s, c, and b. Additional jets are produced in parton 
showering and other parton-level processes considered. Other backgrounds include 
the following: W + jets, Z + jets, £t(172.6 GeV), and WW, WZ, and ZZ vector 
boson production. Both W + jets and Z + jets have exact matrix elements for one 
parton emission and use parton showering for subsequent emissions. 

We would like to apply cuts to these backgrounds without overly diminishing the 
signal (a "good" signal to be quantified in the next subsection on the LHC reach). 
Initially we choose the following cuts labeled CI: 

• n(jets) > 4, 

• $ T > max{100 GeV, 0.2 M eff ), 
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• ErtjiJtJsJt) > 100, 50, 50, 50 GeV, 

• Transverse sphericity: Sx > 0.2, 

where M eff = $ T + E T (ji) + E T (j 2 ) + E T (j 3 ) + E T (j 4 ). After applying the CI cuts 



we arrive at Figures 3.12 and 3.13 Figure 3.12 shows three multiplicity distributions: 



(a) number of jets, (b) number of b-jets, and (c) number of leptons. The kinematic 



distributions in Figure 3.13 are for (a) hardest jet px, (b) second hardest jet px, (c) 
$Xi an d (d) the augmented effective mass, A T (this is the effective mass as earlier, 
but also includes Ex(leptons) and the Ex of all jets). These distributions are used 
to improve on the cuts and each is discussed in the following. 



Figure 3.12 



(a) Jet multiplicity (after relaxing n(jets) > 4) after CI cuts - signals for 
HCAMSBpoints 1 & 2 do not exceed background until very high jet multiplici- 
ties. For instance the HCAMSB1 distribution does not exceed background until 
n(jet) ~ 9. The selection of n(jet) > 4 should be beneficial in this case. 

(b) b — jet multiplicity (after relaxing n(jets) > 4) after CI cuts - signal cross 
section appears harder than background due to the appearance of extra b and 
t quarks from cascade decays as discussed in Section [375 The signals exceed 



background around n(b — jet) ~ 5. We conservatively choose events with at 
least a single b — jet to cut down on background. 

(c) Lepton multiplicity after CI cuts - we can see a much harder signal distribu- 
tion of leptons than for background due to the presence of leptons from cascade 



decays as in the discussion of Section 3.5 The HCAMSB1 signal appears even 



stronger because sleptons are lighter at lower a. Signal exceeds backgrounds 
around n(l) ~ 3, and HCAMSB1 has sufficient strength to be visible here with 
only a few fb" 1 of LHC data [22] |2T]. 
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Figure 3.12: Jet, b-jet, and lepton multiplicity distributions for HCAMSBpoints 1 
(red) & 2 (blue), and SM background (gray), all after CI cuts. 



Figure 3.13 



(a) Et of hardest jet after CI cuts (Et(ji) cut relaxed) - the HCAMSB2 signal 
peaks for ErC/eti) ~ 150 GeV because of b\bi production followed by b\ — > bZ\. 
Signal exceeds background around 450 GeV for HCAMSB1 and 550 GeV for 
HCAMSB2. 



(b) E T of second hardest jet after CI cuts (E T (j 2 ) cut relaxed) - Similar to (a) 
but signal exceeds background at ETijet-i) ~ 350 GeV for HCAMSB1 and ^ 
450 GeV for HCAMSB2. 
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(c) $ T after CI cuts - the HCAMSB2 distribution is softer than for HCAMSB1 
because the $ T originates from third-generation squarks, e.g. bi — > bZ\, whereas 
in the latter case harder Z\s are produced from TeV-scale squarks and gluinos. 

(d) Augmented effective mass - A T = $ T +Y1 ^(jets) + E T {isolated leptons) 
HCAMSB1 has a smooth distribution and exceeds background at At ~ 1600 
GeV. The HCAMSB2 distribution interestingly has two components: a soft 
component with peak at At ~ 750 GeV due to third-generation squark pair 
production, and a hard component with peak at At ~ 2000 GeV due to g and 
qL production. 

And finally, in addition to the distributions just discussed, we consider Figure 



3.14, which shows the oppositely-signed (OS), dilepton invariant mass distribution 
for HCAMSB points 1 & 2 and the background after CI cuts and an extra cut of 
A T > 1500 GeV. This type of distribution is useful in SUSY studies because of the 
appearance of kinematic mass edges due to Z 2 — > or Z 2 — > l^FZi decays. 

For mAMSBmodels the bino-like Z 2 produce a mass edge since it efficiently decays 
to ■ In the HCAMSBmodel however, Z 2 is higgsino-like and rather heavy (see 
Table 3.1 for a comparison). In this case the decay channels are Z 2 — > W^W T , Z t h, 



and ZiZ. Most importantly is the last decay channel which should always be open 
(except when jj, — Y 0, when the a is at its very highest value) and occurs with a 
branching fraction at the tens of percent level. The decay to Z\Z of the bino-like Z 2 
of mAMSBis highly suppressed due to the structure of ZZxZ 2 coupling [20], where 
Z couples to neutralinos only through their higgsino component. Thus we expect 
to be able to distinguish between HCAMSBand mAMSBmodels in the OS dilepton, 
invariant mass distribution because the former has a Z-peak structure and the latter 



does not. This is clear from Figure 3.14, which shows a clear Z-peak structure in 



HCAMSBpoints 1 & 2, and a diminished SM background. 
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Figure 3.13: Kinematic distribution for HCAMSBpoints 1 (red) & 2 (blue) and SM 
background (gray) after CI cuts. 



LHC Reach 

In this section we would like to understand the range of parameters accessible at the 
LHC. To this end, we vary over the range 7713/2 = 30 — 200 TeV and we adopt the 
two following model lines (again, low and high a): 

HCAMSB1: a = 0.025, tan (3 = 10, /1 > 0, 

HCAMSB2: a = 0.15, tan = 10,// > 0. 

We take a somewhat lower value of a in HCAMSB2 than previously because 777.3/2 
could only be extended up to ~ 60 TeV before EW symmetry is improperly broken 
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Figure 3.14: m(l + l~) distribution for HCAMSBpoints 1 (red) & 2 (blue) and SM 
background (gray) after CI cuts. 



see Figure 3.5). The spectrum is shown for the two points above while varying the 



parameter in Figure 3.15 The relative ordering of particle masses are unchanged. 



The distributions of the last section allow us to improve on the background cuts. 
Instead of CI cuts we will now use the following set of cuts that we label C2: 

• n(jets) > 2, 

• Transverse sphericity: St > 0.2, 

• n{b — jets) > 1, 

• E T (j 1 ),E T (j 2 ),$ T > E™\ 

Efj?* will be variable minimum transverse energy value. Parameter space points with 
lower sparticle masses will benefit from a low Etp 1 , while points with high sparticle 
mass but lower production rates benefit from higher E™ 1 . 

We apply the C2 cuts to the following lepton multiplicity channels: Olep, Hep, OS 
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Process 0£ U OS SS U 



QCD(p T : 0.05-0.10 TeV) 
QCD(p T : 0.10-0.20 TeV) 
QCD(p T : 0.20-0.40 TeV) 
QCD(p T : 0.40-1.00 TeV) 
QCD(p T : 1.00-2.40 TeV) 
tt 

W + jets; W -> e,ii,r 
Z + jets; Z — > tt, us 

WW, zz, wz 



73.5 






42.6 


26.5 


37.3 


0.8 


0.6 


0.3 


1253.2 


341.2 


224.9 


60.6 


5.6 


2.8 


61.4 


0.0 


0.77 


0.11 







Summed SM BG 1492.3 374.1 266.1 0.26 0.25 



HCAMSB1 
HCAMSB2 



100.1 
223.5 



53.2 
58.7 



13.1 
4.6 



2.4 
1.7 



3.3 
0.35 



Table 3.4: Estimated SM background cross sections (plus two HCAMSB benchmark 
points) in fb for various multi-lepton plus jets +$t topologies after cuts C2 with 



E C T = 100 GeV. 



dilep, SS dilepQ 31ep, and 41ep. As a supplement to the above cuts, in the Olep channel, 
it is required that the .^-(nearest jet E?) transverse opening angle is constrained to 
30° < Acj)($ T , E T (j nearest )) < 90°. Furthermore, all isolated leptons are required to 
have a minimum of 20 GeV, and events with single leptons are required to have 
transverse mass M T {l,$ T ) > 100GeV to reject leptons from W decays. 

The results of C2 cuts on SM background and the two model lines is shown in 



Table 3.4 for E^ = 100 and for 2M events. Backgrounds cannot be detected for several 

multi-lepton channels. We consider an observable signal to be one that satisfies the 
4 We do not require "same flavor" in the dilepton channels. 
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following criteria for an assumed integrated luminosity: 

i. S/B > 0.2, 

ii. Signal has at least 5 events, and 
Hi. S > b\f~B, 

Where S and B are the number of events for signal and background respectively. The 
first requirement is imposed to prevent the use of data that has a small signal on top 
of a large background, and that requires the use of extremely precise backgrounds 
[23] . The last requirement means that we require the signal to have a 5cr statistical 
significance above the background. 

Using the criteria above and assuming 100 fbT 1 of data, we plot signal and back- 
ground versus m 3 / 2 for C2 cuts with each of Ej, = 100, 300, and 500 GeV and in each 



of the lepton multiplicity channels. Figure 3.16 shows the 5<r/5 event SM background 



values with dotted horizontal lines, while the signal is represented by solid lines. Sig- 
nal values above the horizontal lines qualify as observable. The value of 7713/2 where 
the solid line meets its corresponding (same color) dotted line is the value of m 3 / 2 
that is accessible at the LHC. A SUSY signal is not observable where the solid line 
appears below the horizontal line. 



The summary of the results from the plots in Figure 3.16 is given in Table 3.5 
In the table, the upper entries are the m 3 / 2 reach values for point 1 (a = 0.025) and 
the lower entries are the same for point 2 (a = 0.15), and each row corresponds to 
a different Ej, value. The highest gravitino mass probed for Ej, = 100 GeV is 80 
TeV and is found in the 31 channel. However, we can see by going to the harder 
cut, E^p = 500 GeV, in the 01 channel, the reach is extended up to 115 TeV for m 3 / 2 
for HCAMSB1. This value of m 3 / 2 corresponds to a gluino mass of nig ~ 2.4 TeV. 
Similarly, the best reach for HCAMSB2 is in the 0/ channel with Ej* = 500 GeV, with 
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Table 3.5: Estimated reach of 100 fb" 1 LHC for m 3/2 (TeV) in two HCAMSB model 
lines: a = 0.025 (upper entry) and a = 0.15 (lower entry), in various signal channels. 

the 105 TeV gravitino mass being probed, corresponding to a gluino mass of ~ 2.2 
TeV. 

3.7 Summary 

In this chapter we have given a complete analysis of the HCAMSBmodel relevant to 
physics at the LHC. We calculated the spectrum from which we were able to further 
calculate branching fractions and production cross sections, and to simulate LHC 
collider events for comparison again SM backgrounds. The final result of our analysis 
is that after 100 fb _1 of data collection running at 14 TeV, the gravitino mass reach is 
m 3/2 ~ H5 TeV {nig ~ 2.4 TeV) for low values of a, and m 3 / 2 ~ 105 TeV {nig ~ 2.2 
TeV) for large a. Since the spectra is so similar for /i > and fi < 0, we do not expect 
changes in the reach due to the change in the sgn{u). If instead tan /3 was moved 
to higher values, we do not expect there to be differences in the 0/ and 11 channels 
where we have maximal reach. Although it is possible there would be differences in 
the multi-lepton channels due to enhanced -ino decays to rs and 6s at large tan j3. 
In references [16j[30j, the mAMSBreach was shown to be around m g ~ 2.75 TeV for 
low values of m , which is somewhat larger than for HCAMSB. This is due to the 
fact that squark masses tend to cluster around a common mass scale mo in mAMSB, 
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C2 Cuts + 01 C2 Cuts + 11 




C2 Cuts + 31 C2 Cuts + 21 (OS) 




m rTeV] m rTeV] m [TeV] m,_ [TeV] 




Figure 3.16: Cross section versus m 3 / 2 for HCAMSBpoints 1 & 2, after C2 cuts and 
in various lepton channels. 
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while in the HCAMSBcase right and left squark states are highly split. 

The mAMSBand HCAMSBmodels appear to have similar signatures at the LHC. 
Both would have multi-jet, multi-lepton, and $ T events. In addition to these, long 
chargino tracks of 0(1 — 10) cm are possible. However, there are major differences 
between these two models that are summarized below. 

• HCAMSBhas a highly left-right split scalar spectrum due to the left- right 
asymmetry in hypercharge assignments. For mAMSBthere is instead a near- 
equality between left and right mass parameters. The lightest stau, fx, is mainly 
left-handed, while in mAMSBit is mainly right-handed. While it is possible that 
the left-right mixing could be determined at the LHC, this task would be easily 
done at a linear collider with polarized beams |15j . 

• We saw that high rate production of 6s and is was due to light t\ and b\ states 
relative to the gluino. In particular, with the HCAMSBmodel we expect high 
multiplicities of final state b-quarks, whereas for mAMSBmodels much few are 
expected (this also depends on the chosen tan (3). 

• We also see that the ordering of \i > Mx > M 2 implied a bino-like Z 2 for 
mAMSB, while for HCAMSB, with M x > fi > M 2 , it is higgsino-like. This 
crucially changes the decays of Z 2 and leads to an important distinction between 
mAMSBand HCAMSBmodels: the m(ll) distribution for OS, dilepton events 
has a Z — > peak in the latter, while the former would have a smooth 
distribution with no peak (except at large tan where there is greater mixing 
in the neutralino sector). Thus, HITs (see Chapter [2j from charginos and a 
Z-peak from Z 2 decays could yield a promising signal for HCAMSBand LHC. 
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4 



inoAMSB 



4.1 Introduction to inoAMSB Models 



The combination of anomaly mediated (Section lL4) and gaugino mediated [67J soft 



SUSY breaking parameters is suggested by two classes of models that arise in type 
IIB string theories compactified on a Calabi-Yao Orientifold (CYO). These model 
types differ in the number of Kahler moduli and they are: 

• Single Kahler modulus (SKM models) - KKLT type with uplift of soft terms 
and cosmological constant coming from one-loop effects, and generally favors 
smaller CY compactification volumes; 

• Large Volume Scenario (LVS models) - require at least two moduli and, as the 
name implies, favors larger CY compactification volumes. 

The F-terms of moduli are responsible for SUSY breaking. In both cases the moduli 
are stabilized using a combination of fluxes and non-perturbative effects while the 
interactions with the MSSM are gravitational. The models also share two crucial 
features: gauginos receive their masses only from Weyl anomalies (string contributions 
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are suppressed) and the soft scalar and trilinear parameters are naturally suppressed 
as usual in gaugino mediated scenarios (see [HZ] [20] ) • 

Having just analyzed the H CAM SB mo del, it should be possible to study other 
models in a similar fashion. In this chapter we do just that for a class of string 
models whose soft parameters arise in a pattern that we call the "Gaugino Anomaly 
Mediated SUSY Breaking" pattern, or more succinctly, "inoAMSB" [26J. The string 
theory origin for this model is quite different than for the case of HCAMSB, but we 
will see some similarities in the spectra that will lead to collider signatures already 
familiar from the last chapter. Furthermore, we will more or less follow the same 
procedure for arriving at the collider analysis level. Schematically this procedure is 
as follows: 

• define the model high-scale boundary conditions, 

• discuss the evolution of theory parameters from high-scale to TeV scale, 

• explore the allowed parameter space subject to experimental constraints, 

• choose the generally representative points for collider analysis and examine 
branching fractions and cross sections to determine possible signatures, 

• and finally, determine the parameter reach for an experiment to produce a 
statistically significant signal. 

4.2 Setup of the inoAMSBModels 

Before moving on to the details of the models we should first consider a few prelimi- 
naries. A supergravity model is defined in terms of three functions: a superpotential, 
a Kahler potential, and the gauge kinetic function. They each have implicit $ depen- 
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dence and can be expanded in a series of MSSMsuperfields (C a ) as follows: 

W = W($) + n($)H A H u + l -Y aM C a CP& + ... (4.1) 
K = K($, $) + K a p(<5>, $) C a C? + [Z(<t>, $)if d if u + /i.e.] + . . . , (4.2) 

fa = /«(*). (4.3) 

The role of these moduli fields is to break supersymmetry. This happens once they 
are stabilized and acquire VEVs such that at least one of them has a non-zero F-term. 

Each of the first terms in the above expansions have implicit dependence on the 
Calabi-Yau Orientifold (CYO): 

K = -21n(V +\)- ln(i j H A Q(U, U)) - ln(S + S), (4.4) 

W = f G 3 An + J2 A ae' atTl - (4.5) 
J i 

Most notably, the first term in K contains the CYO volume, V, and a stringy cor- 
rection £ that depends on the Euler character of the CYO and the real part of the 
dilaton superfield S. The last term in W is non-perturbative and is responsible for 
stabilizing the Kahler modulus at the minimum of the potential as in KKLT models. 
Q are 3-forms on the CYO and are functions of U r , the complex structure moduli. 
The standard SUGRA potential is given by 

V($) = F A F S K AB - 3|m 3/2 ($ )| 2 , (4.6) 

where the SUSY breaking F-terms at the minimum are given by F A = e K l 2 K AB D^W , 
and the covariant derivatives are given by Da = <9a + Ka, where subscripts indicate 
derivatives with respect to fields. Also, $ stands collectively for all moduli that de- 
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scribe the internal geometry of the CYO plus the axion-dilaton. The minimum of the 
potential, Vo, gives the cosmological constant (CC) of the theory. The CC is known 
to be small, and taking it to be close to zero implies that at the minimum of the 
potential F A ~ m 3 / 2 . 



Single Kahler Modulus Scenario (SKM) 

This model has IIB string theory compactified on a CYO, has only one Kahler mod- 
ulus, T, and the MSSMfields come from open string fluctuations on a stack of D3 
branes. Other moduli and the axio-dilaton are stabilized by internal fluxes and non- 
perturbative effects, as in KKLT . 

At the classical level the CC would be small in magnitude but negative and the 
soft masses would be highly suppressed. However, because SUGRA theories have 
quadratic divergences at the quantum level (unlike broken global SUSY ), there are 
corrections to these quantities that are dependent on the string scale cutoff A. Such 
contributions can serve to uplift the CC to small and positive values and can gener- 
ate soft SUSY breaking masses proportional to ^7713/2- The cutoff is taken to be the 
string scale and can be between ~ 10 14 GeV and ~ M GUT . 

Large Volume Scenario (LVS) 

In this class of models |28j, IIB string theory is again compactified on a CYO. However, 
now more than one Kahler modulus is considered, Tj (i — 1, . . . , hu). In particular, 
in the simplest situation there is a large modulus, r^, and small moduli, (r s , r a ), con- 
trolling the overall size of the CYO and the volume of two small 4-cycles respectively. 
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The total volume is then given by 

V = r b 3/2 - ~ rT- (4.7) 

This is referred to as a "Swiss Cheese" model. Again the MSSM may be located on 
D3 branes at a singularity. Alternatively, it could be placed on a stack of D7 branes 
wrapping a four cycle (taken to be the one labelled by the index "a"). In this case, 
it has been argued [35] that the necessity of having chiral fermions on this brane 
prevents this cycle from being stabilized by non-perturbative effects and it shrinks 
below the string scale. Effectively, this means that the physics is the same as in the 
D3 brane case. 

Extremizing the potential leads to an exponentially large volume [2S] V ~ e aTs , r s ~ 
£. It turns out that the suppression of FCNC effects lead to V ~ 10 5 /p [U] (where 
lp is the Planck length), so the string scale is M string ~ Mp/y/V ~ 10 155 GeV. The 
minimum of the potential (CC) is given by Vq ~ ~~ i n ^ 2 y - This minimum can be 
uplifted to zero when S and U r acquire (squared) F-terms of the order ln 3 ^ 3/2 v • Clas- 
sical contributions to the scalar and slepton masses are also of this same order. With 
the above lower bound on the volume, this means that even for m 3 / 2 ~ 100 TeV, the 
classical soft terms are ~ 100 GeV. Of course if one wants to avoid fine-tuning of the 
flux superpotential, it would be necessary to take even larger values of V correspond- 
ing to a string scale of 10 12 GeV. In this case the classical soft terms are completely 
negligible (for m 3 /2 ~ 100 TeV) but the (classical) /x-term is also strongly suppressed. 

In the rest of this section the holomorphic variable associated with the large mod- 
ulus u will be called T. 
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Gaugino Masses - Weyl Anomaly Effects 

For a generic version of supergravity, the gaugino masses satisfy the following relation 
at the classical level: 

M. 1 

F A d A f a {®). (U 



9l 2 



In the single Kahler modulus model the MSSM resides on D3 branes at a singularity. 
In the LVS case, we may either have the MSSM on D3 branes at a singularity, or we 
may have it on a stack of D7 branes wrapping a four cycle which shrinks below the 
string scale. In both cases the classical gauge coupling function is effectively of the 
form 

fa = S (4.9) 

where S is the axio-dilaton (for more details see (47] section 3.1 and references therein). 
The important point is that it is independent of the modulus T. So at a classical 
minimum, where the SUSY breaking is expected to be in the T modulus direction, 
the string theoretic contribution to the gaugino mass is highly suppressed. 

However, there is an additional contribution to the gaugino mass due to the (super) 



Weyl anomaly as discussed in Section 1.4 This comes from the expression for the 



effective gauge coupling superfield that has been derived by Kaplunovsky and Louis 
[HE] (KL f\ For the gaugino masses, the relevant contribution comes from taking the 
F-term of 

3c ai _ T a (r 



H a ($,T,r z ) = / a ($) - — |lnC - -i^r z . (1.37) 



i7T 



4vr 2 



Here, the first term on the RHS is the classical term; the second comes from the 
anomaly associated with rotating to the Einstein-Kahler frame. c a = T(G a ) — 
^ r T a (r) is the anomaly coefficient and the last term comes from the anomaly asso- 
ciated with the transformation to canonical kinetic terms for the MSSM fields. Also 



1 As explained in [46], the usual formulae for AMSB need modification in the light of [68 . 
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note that we have ignored the gauge kinetic term normalization anomaly P, HH] which 
is a higher order effect. The chiral superfields 0, r that generate these transforma- 
tions are given by, 

lnC + lnC = -K| harmonic , (4.10) 
r z + f z = IndetAg. (4.11) 

The instruction on the RHS of the first equation is to take the sum of the chiral and 
anti-chiral (i.e. harmonic) part of the expression. After projecting the appropriate F 
terms we arrive at the following expression: 



^ = F*ajJ. - £F*K A - £ ^F*d A Indet K^. (4.12) 



The first (classical) term is greatly suppressed relative to m 3 / 2 since the T-modulus 
does not contribute to the classical gauge coupling function as discussed earlier (see 



paragraph after (4.8)). The dominant contribution therefore comes from the last two 
(Weyl anomaly) contributions. It turns out that (after using the formulae F T = 
— (T + T) 7713/2, Kt = — 3/(T + T) and K a p = k a p/(T + T) which are valid up to 
volume suppressed corrections), this yields^J 



M * = Tirl™3 /2 , (4.13) 

1D7T 



where b a = — 3T(G a ) + J2 r T a (r) is the beta function coefficient. 



2 Note that we expect the Weyl anomaly expressions for the gaugino masses given below to be valid 
only because of the particular (extended no-scale) features of this class of string theory models. It so 
happens that these are exactly the same as the expressions given in what is usually called AMSB: but 
that is an accident due entirely to the fact that in these extended no-scale models the relationship 
F a Ka — 3m 3 / 2 is true. 
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Scalar Masses, Trilinear Couplings, /i and Bfi terms 

Here are summarized the results from this class of string theory models for the values 
of the soft parameters at the UV scale, i.e. A ~ M string ~ Mp/y/V. These values 
should be the initial conditions for the RG evolution of these parameters. In the LVS 
case, it was estimated [UJ that the lower bound on the CYO volume was V > 10 5 . 
Also, typical values of h 2 i ~ O(10 2 ) are chosen for the number of complex structure 
moduli. The gravitino mass is chosen to be m 3 / 2 ~ |W|Mp/V ~ 50 TeV. Such a 
large value of m 3 / 2 allows us to avoid the SUGRA gravitino problem, which leads 
to a disruption of Big Bang nucleosynthesis if m 3 / 2 ~ 5 TeV and Tr ~ 10 5 GeV 

UJHUIHI]. 

Unlike the gaugino masses, scalar masses and trilinear soft terms do not acquire 
corrections from the Weyl anomaly. They are essentially given at the UV scale by their 
classical string theory value plus one loop string/ effective field theory corrections. In 
the hn = 1 case, the classical soft terms are essentially zero while in the LVS case 

m o ~ 1 , i m3/2 ) , — ~\A2i% A) < rriQ. (4.14) 

y yin m 3/2 V J n 

As discussed in [UJ the estimate of the fi, Bfi term comes from the generic case 
of assuming that the uplift of the negative CC is distributed amongst the complex 
structure moduli as well as the dilaton. It was estimated in that reference that the 



F-term for a complex structure modulus is \F m \ ~ m 3 / 2 / \fh2~1 (In m 3/2)V- Since the 
expression for the /1 term involves a sum over h 2 i terms this gives the estimate in 



(4.14). The parametric dependence on the volume factor in the upper bounds for the 
fi, Buffi terms goes as \/h 2 \l W and clearly favors the "small" values CYO volume 
i.e. V ~ 10 5 and the large values of the number of complex structure moduli (~ 10 2 ). 
However it should be stressed that these estimates are on a different footing than 
the rigorous calculation of m and the gaugino masses, and indeed it is possible that 
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these LVS models may have a /i-problem. 

After adding quantum corrections at the UV scale, both cases give similar values 
for the soft terms. As an example, this is illustrated for two values for the CYO 
volume: 

• V ~ 10 5 , M string ~ A ~ lQ- 2 - 5 M P ~ 10 15 - 5 GeV. Then, 



H ^ < 250 GeV, m ~ 25 GeV, A < m . (4.15) 



V ~ 10 12 , M string ~ A ~ 10~ 6 M P ~ 10 12 GeV. Then, 



H ^ EH < io- 1 GeV, m ~ 10~ 2 GeV, A < m . (4.16) 



The second very large volume case can be accessed only in the LVS model. Our fa- 
vored case however is the first one with a volume around V ~ 10 5 in Planck units. It 
should be noted that even with this value, it is only expected that stringy corrections 
to the numerical estimates of the above values of the soft masses etc. will be of order 

a'/R 2 CYO ~ VV 1/3 ~ 10%. 

The first case is at the lower bound for the volume. This gives the largest allow- 
able string scale. This is still somewhat below the apparent unification scale, but it 
is close enough that (allowing for undetermined 0(1) factors) the GUT scale may be 
used as the point at which to impose the boundary conditions. This is useful for the 
purpose of comparing with other models of SUSY mediation where it is conventional 
to use the GUT scale. 

The second case above corresponds to choosing generic values of the flux superpo- 
tential, while the first needs a fine tuned set of fluxes to get \W\ ~ 10~ 8 , in order to 
have a gravitino mass of ~ 10 2 TeV, though in type IIB string theory general argu- 
ments show that there exist a large number of solutions which allow this. The most 
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significant problem with the second case (apart from the fact that there is no hope of 
getting a GUT scenario) is the extremely low upper bound on the \x term. In other 
words, there is a serious \i- problem. The first case also may have a \i term problem, 
but again since these estimates are accurate only to 0(1) numbers, it is possible to 
envisage that the problem can be resolved within the context of this model. 

In any case, the string theory input is used to suggest a class of phenomenological 
models. Given that in both the GUT scale model and the intermediate scale model, 
the soft scalar mass and A term are suppressed well below the weak scale and assumed 
to be zero for these at the UV scale, while the gaugino masses at this scale are given 



as was the case in the last subsection. 

The case when the input scalar mass m is non-negligible will also be discussed. 
This would be the case for instance in the SKM model with smaller volumes and/or 
larger values of /121, and also in the case of LVS with the volume at the lower bound 
but with larger values of h 2 i- 



4.3 Spectrum, Parameter Space, and Constraints 



by 



M a = 




(4.17) 



The parameter space for this model is 



m 3 / 2 ,tan j3, and sgn(/i), 



and, given the discussion of the previous section the boundary conditions 



m = A 



— 







(4.18) 
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are imposed along with AMSB masses for gauginos at M str ing/GUT- Note that even 
when m and A are negligible at the GUT, scalars will be massive at the weak scale 
due to uplift from gauginos in the renormalization group running. 

To illustrate the running of soft masses from the GUT scale to the weak scale, 
Figure 4.1 shows gaugino masses (top) and scalars (bottom) for the point m 3 / 2 = 50 
TeV, tan (3 = 10, and \i > for M str i ng = Mqut- Gauginos begin with the usual ratio 
of (Mi : M 2 : M 3 ) = (6.6 : 1 : -4.5), but at the weak scale the 1 and 3 masses switch 
in ordering leaving |Ms| >> Mi > M 2 . The lower value of M 2 indicates that the 
lightest neutralino will be wino-like. Scalars begin with negligible masses as required 
by the boundary conditions. However the mass parameters are uplifted by radiative 
self-couplings. In particular, the large value of M\ is responsible for the uplift near 
the Mqut- In the early running the right-slepton, E 3 , moves to the highest values 
because it has the largest hypercharge (Y=2). Because the hypercharges for the left- 
sleptons are lower than for the right-sleptons, they appear to have lower masses at the 
weak scale, which leaves us with a left-right split spectrum for sleptons. This effect 
is similar to what was seen in the case of HCAMSB, although the splitting for that 
model was more severe because of the much heavier M\. Squarks receive extra uplift 
from QCD effects leaving them in the TeV range, while sleptons are generally lower: 



approximately between 200-400 GeV. Also, since Z\ is wino-like and M 2 < w m | 1 
charged LSPs are not of concern. 

At the weak scale, physical masses and mixings are computed as usual and 1-loop 
corrections are added. The resulting spectrum for two inoAMSBpoints with m 3 / 2 = 
50 and 100 TeV is listed in columns 3 and 4 of Table 14.11 In column 1 of the same 
table is the mAMSBpoint with tuq = 300 GeV in column 2 is an HCAMSBpoint with 
a = 0.025 for comparison. By inspection it is seen that the features of the previous 
paragraph are represented in the table: large left-right splitting of sleptons, generally 
heavier squarks, and a Z\ with a mass close to M 2 reflecting that it is wino-like. By 
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Q (GeV) 



Figure 4.1: inoAMSB soft SUSY breaking parameters as a function of energy scale Q 
for model parameters m 3/ / 2 =50 TeV, tan (3=10, /i > 0, and M str - ing =MQ UT . 



comparison, HCAMSBhas a much larger left-right splitting in the sleptons (reflective 
of the heavier M\ in that case), whereas mAMSBhas a much smaller splitting (nearly- 
degenerate with mg R < rrie L ). 

The inoAMSBparameter space is simpler than either the mAMSBor HCAMSBcase 
because it has one less model parameter. We begin by recognizing that once a point is 
chosen, the increasing 717,3/2 value increases all masses but leaves the relative hierarchy 



unchanged. This is demonstrated in Figure 4.2 (a) where the spectrum is plotted 
versus m 3 / 2 for tan (3 = 10 and > 0. Everywhere in the parameter space of 
the inoAMSBmodel the hierarchy is generally in the order rrig > rriq > > 
> m e L ,u L > m z 1 ,w 1 - We then expect squark pairs to be produced in LHC events 
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Table 4.1: Masses and parameters in GeV units for four case study points mAMSB 1, 
HCAMSB1, inoAMSBl and inoAMSB2 using Isajet 7.80 with m t = 172.6 GeV and 
> 0. Also listed are the total tree level sparticle production cross section in fb at 
the LHC. 



either directly or through the decay of pair-produced gluinos. Furthermore, because 
squarks are heavier than sleptons we can expect squark cascade decays to two hard 
jets + isolated leptons through q — > qZi — > q^l T . 
In frame (b) of Figure 



4.2 



taking m 3 / 2 = 50 TeV, an example of the spectrum 



versus tan fi is given. Some interesting effects occur in the region of high tan fi. There 
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are noticeable dips in the masses for b and f eigenstates since larger tan (3 results in 
larger Yukawa couplings for these particles, and thus greater suppression of masses 
when approaching the TeV scale. The tua parameter decreases most significantly in 
the plot. Because of the increased down-type Yukawa couplings for large tan (3, mn d is 
pushed towards negative values. Since the EW-breaking minimization conditions give 
m\ ~ m# ~ m2 H u i the pseudoscalar too decreases for larger tan (3. Negative values 
of rn 2 H d signal improper breaking and this occurs around tan (3 ~ 42. And finally, 
rriA cannot fall below the LEP2 bounds on the Higgs. In the plot rrih ~ 111 GeV 
which we consider near the acceptable edge due to the ±3 GeV error in the theory 
calculation. As tan (3 increases, the LEP2 bound is violated, and shortly thereafter 
the EW symmetry is not properly broken. It can also been seen that at the lowest 
tan (3 values, the lightest Higgs is too light. 



With the features of Figure 4.2 in mind, we look at the entire tan (3 — m 3 / 2 



parameter space in Figure 4^3 The orange region is specifically excluded by LEP2 
bounds on the chargino mass, and all m 3 / 2 values between this region and zero are 



likewise excluded. The brown regions are also excluded as in Figure 4.2 (b): the high 
tan (3 regions has improper EW-breaking and the lower region (just above tan (3 = 0) 
is where the light Higgs is too light. Gluino contours are also shown in the plane, 
and it is seen that the major (quasi-linear) dependence is on m 3 / 2 - A gluino with the 
higher masses shown will not be accessible at the LHC. 

Until now we have made the assumption that M str i n g = Mqut, but this is not 
necessarily so as pointed out earlier. In fact, we should expect that M str i n g will be 
somewhat lower than the GUT scale in order to have enough FCNC suppression. 



Figure AA is a plot of the spectrum as a function of the string scale, and the latter is 
taken as low as 10 11 GeV. As M string is lowered the spectrum is seen to spread out, 
and there are cases of mass re-ordering. The most important feature of this plot is 
that around 5 x 10 13 GeV, the LSP is the tau-sneutrino. There are severe limits on 
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Figure 4.2: Sparticle masses: (a) as function of m.3/2, with tan (3 = 10, fi > 0, 
and Muring = Mqut] (b) as a function of tan (3, with m 3 / 2 = 50 TeV, \i > 0, and 

^string = M GUT . 



stable sneutrino dark matter [51. Then, in order to avoid these limits and to have a 



Zi LSP, we do not consider such low M. 



string- 
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Finally, we also have been assuming that tuq = 0. However this is not completely 
realistic, as they are only expected to be suppressed. The spectrum is shown again 



in Figure 4J3 this time as a function of mo, the classical mass added universally to 
all scalars. We adopt values 777,3/2 = 50 TeV and tan (3 = 10 for this plot. As m 
increases beyond zero, it is seen that the spectra change little so long as m ~ 100 
GeV, and also the mass orderings remain intact. For larger values of mo, the left- 
and right-slepton masses begin to increase, with first surpassing m^ 2 , and later 
even surpasses m^ 2 . At these high values of m , decay modes such as Z 2 — > 
IF would become kinematically closed, thus greatly altering the collider signatures. 
However, generically in this class of models, we would not expect such large additional 
contributions to scalar masses. 




Figure 4.3: tan (3 — m 3 / 2 parameter space, with /j, > and M str i ng = Mqvt- Also 
shown are gluino mass contours in the parameter plane. 
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Figure 4.4: inoAMSBspectrum versus string scale for m 3 / 2 = 50 TeV, tan (3 = 10, 
and n > 0. 
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Figure 4.5: inoAMSBspectrum versus classical common soft mass, rriQ. The parame- 
ters for this plot are m 3 / 2 = 50 TeV, tan (3 = 10, and fi > 0. 
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BF(b s 7 ; and (g-2)^ 

Along with experimental constraints on the inoAMSB models from LEP2 limits on 
and rriyy , there also exist indirect limits on model parameter space from comparing 
measured values of BF(b — > 57) and Aa M = (g — 2)^/2 against SUSY model predic- 



tions. Figure 4.6 (a) shows regions of the branching fraction for BF(b — > 57) in the 
inoAMSB model versus m 3 / 2 and tan (3 variation, calculated using the Isatools sub- 
routine ISABSG[13j. The red-shaded region corresponds to branching fraction values 
within the SM theoretically predicted region BF{b — > sj)sm — (3.15 ± 0.23) x 10~ 4 , 
by a recent evaluation by Misiak|84J). The blue-shaded region corresponds to branch- 
ing fraction values within the experimentally allowed region [29]: here, the branching 
fraction BF(b — > 57) has been measured by the CLEO, Belle and BABAR collabo- 
rations; a combined analysis (29] finds the branching fraction to be BF(b — > 57) = 
(3.55 ± 0.26) x 10 -4 . The gray shaded region gives too large a value of BF(b — > 57). 
This region occurs for low m 3 / 2 , where rather light t\ and W\ lead to large branching 
fractions, or large tan j3, where also the SUSY loop contributions are enhanced 



Figure 4.6 (b) shows the plot of the SUSY contribution to Aa M : Aa^ USY (using 



ISAAMU from Isatools [17J). The contribution is large when 777,3/2 is small; in this 
case, rather light fx^ and masses lead to large deviations from the SM prediction. 
The SUSY contributions to Aaf, USY also increase with tan B. It is well-known that 
there is a discrepancy between the SM predictions for Aa M , where r decay data, used 
to estimate the hadronic vacuum polarization contribution to Aa M , gives rough ac- 
cord with the SM, while use of e + e~ — > hadrons data at very low energy leads to a 
roughly 3<r discrepancy. The measured Aa M anomaly, given as (4.3 ±1.6) x 10~ 9 by 
the Muon g — 2 Collaboration [39], is shown by the black dotted region. 
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Figure 4.6: (a) BF(b ->■ 57) and (b) SUSY contribution to (g - 2)^ in the inoAMS- 
Bparameter plane with M string = Mqut- 
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mAMSB HCAMSB1 inoAMSBl inoAMSB2 



a[fb] 


7.7 x 10 3 


7.4 x 10 3 


7.5 x 10 3 


439 


g, q pairs 


15.0% 


15.5% 


19.1% 


3% 


EW — ino pairs 


79.7% 


81.9% 


75.6% 


93% 


slep. pairs 


3.7% 


0.8% 


3.1% 


3% 




0.4% 


0.2% 


0.1% 


0% 



Table 4.2: Cross sections for pp — > SUSY at 14 TeV. 



4.4 inoAMSBat the LHC 



Sparticle Production 



Table 4^ shows the distributions in the cross section of various sparticle pair-production 
channels. EW-ino pairs dominate over all other forms of production. As in the case 
of HCAMSB, because Z\ is stable and W x decays to soft ^ and Z\, the processes 
pp — > W^W^ 1 , WfZ\ do not produce sufficient visible energy to meet detector trigger 
requirements. We must instead rely on detecting produced squark and gluino pairs. 
We will will be mainly interested in pp — > gg, qg, qq because of their significant 



cross sections and visible final states. Figure 4.2 shows that there are significant cross 
sections for slepton production as well, but those cross sections are still lower than 
QCD pairs and the LHC reach in 1713/2 is much less. Therefore, for this analysis we 
focus strictly on QCD pair-production. 

Since sparticle masses depend mainly on m 3 / 2 , the plot of multiple pair-production 
cross sections is shown as a function of this parameter in Figure |4.7 We see that for 



m 3/2 ~ 65 TeV gq and qq rates are comparable, but for higher m 3 / 2 that level squark 
pairs dominate. Rates for gg can be significant, but they are never dominant. 
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Figure 4.7: Production cross sections for various pairs versus for tan (3 = 10, 
and fi > 0. 

Sparticle Decay 

We saw that the inoAMSBRGE running leads to a hierarchy of masses in which all 
squark masses are close in value and less than the gluino mass. Thus, produced 
gluinos decay nearly equally to squark species through g — > qq. Left squarks decay 
mainly through q^ — > q + wino, and the calculation shows q^ — > q'W\ at 67% and 
q~L qZ\ at 33% over the entire parameter space. Right squarks on the other hand 
decay to q + bino, and therefore the decays qn —> qZ 2 occur at 97% over all the 
parameter values. 

The RGE evolution also led to a large left-right splitting in the sleptons, and this 
should have noticeable effects. Left sleptons decay mainly to wino + lepton and the 
calculation shows II — > IZ\ at 33% and II —> v\iW\ at 67% over all parameters. 
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Sneutrinos decay invisibly to v\Z\ at 33%, but has 66% visible branching to IW\. If 
right-sleptons were heavy enough, they would decay to bino + lepton, but since they 
are too light they decay via Ir — > eZ\ at 78%, as well as through three-body modes 
l R — > /~T + ff and l~r~fi at 13% and 7% respectively. 



inoAMSBLHC Events 

Isajet [SB] was used for the LHC event simulation, the details of which can be found 



in Section L5 Two-million events for inoAMSBpoint 1 of Table 44 were generated in 
addition to the same Standard Model backgrounds discussed in the Chapter |3j which 
are W + jets, Z + jetfl ti, WW, WZ, and ZZ. 

We first apply a rudimentary set of cuts on signal and background that we will 
label CI. This first round of cuts will help to understand the underlying properties 
of LHC events and will give us information on how to make better cuts. The CI cuts 
are as follows: 

• n(jets) > 2 

• $ T > max(100 GeV,0.2M e// ) 

• E T ( 3u j 2 ) > 100,50 GeV 

• transverse sphericity St > 0.2, 

where M eff = $ T + £ T (ji) + E T (j 2 ) + E T (j 3 ) + E T (j 4 ). 

Since the sparticle production is dominated by gq and qq we expect to have two 
very hard jets in each event through the squark decays to EW-ino + quark. This 



is seen in Figure 4^ where the first and second hardest px distributions are seen to 
emerge from background around 450 GeV and 250-300 GeV respectively. 

■ , -\jets > isolated 

Figure 4.9 shows the $ T and At(= $t + 2^ &t + Yl leptons Et) distributions 



3 Vector boson + jets uses exact matrix element for a single parton emission and parton showering 
for subsequent emissions. 
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Figure 4.8: pt distributions for hardest (left) and second hardest (right) jets for 
inoAMSBl and SM backgrounds after CI cuts for 14 TeV pp collisions. 



for inoAMSBl with SM backgrounds. The $ T distribution emerges from background 
around 500 GeV and appears rather hard because squarks have a significant invisible 
branching to Z\&. The At distribution has a peak around 400 GeV that is buried 
under background, and a broad hard peak at higher At- The softer peak originates 
in chargino, neutralino, and slepton pair production, while the hard peak is due to 
squark and gluino pair production and emerges from background around 1400 GeV. 



We also examine jet (left) and lepton (right) multiplicities in Figure 4.10 We expect 
to see hard dijet events from squark decays, but the jet distribution shows a broad 
peak between rij = 2-5 as well as much higher values. These extra jets occur because 
of cascade decays and ISR and because in this first round of cuts $ T > 50 GeV which 



is rather low compared to when the hardest jets exceed background in Figure 4.8 



The leptons are seen to stand out above background already at 3 leptons, and so we 
see already that with the minimal cuts a signal with 3 leptons, at least 2 jets, and 
$ T should appear well above background. 

Again, as was the case for mAMSBand HCAMSBalready encountered, HITs (see 
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Figure 4.9: Et and At distributions for inoAMSBl and SM backgrounds after CI 
cuts for 14 TeV pp collisions. 

Section [2]) from long-lived charginos is an important signal for models with wino-like 
Zis, particularly AMSB models. If the anomaly mediation supersymmetry breaking 
pattern is the correct one, we expect HITs to play an important role in the discovery 
of SUSY at the LHC. After an AMSB discovery, deciphering further whether a model 
is correct would be the next step. After discussing the reach in the next subsection, 
it will be shown that if mAMSB, HCAMSB, and inoAMSBare accessible at the LHC, 
it would be possible to unambiguously distinguish between all three models! This is 
a major result of this work. 



LHC Reach 

Now we would like to find what are the maximum inoAMSBparameter values ac- 
cessible at the LHC after one year of data collection (100 fb" 1 ) at y/s = 14 TeV. 
This amounts to finding the reach in the parameter because it makes the most 



significant contribution to sparticle mass as seen in Figure 4.2 
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Figure 4.10: Jet multiplicity (left) and lepton multiplicity (right) distributions for 
inoAMSBl and SM backgrounds after CI cuts for 14 TeV pp collisions. 



Taking what was learned in the last section by the CI cuts on signal and back- 
ground, we adopt the following new set of cuts labeled C2: 

• n(jets) > 2 

• S T > 0.2 

where Ej, is variable in order to maximize the reach. Parameter space points with 
lower masses benefit from lower E^. Heavier particles have lower cross sections but 
also have higher energy release per event and thus are more visible for higher Ej,. 
Additionally, we apply extra cuts for multilepton channels. For the 0/ channel we 
apply a cut on the transverse opening angle between $ T and the nearest jet, 30° < 
A<f)($ T ,E T (j near ) < 90°). For all isolated leptons a minimum of 20 GeV is re- 
quired. For the sake of brevity we will consider 0/, 21(0 S), 31, and 4/ channels because 
they provide the best reach. 

After applying C2 cuts with E^ = 100 GeV on SM backgrounds (2M events each) 
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22.6 


8.4 


1.3 




Z + jets; Z —)• rf , z/s 


752.9 


11.1 


1.3 


0.2 




WW, zz, wz 


3.4 


0.3 


0.25 







Summed SM BG 


4776.1 


1697.8 


467.7 


181.9 


9.4 


inoAMSBl 


112.7 


85.7 


27.6 


36.0 


7.5 



Table 4.3: Estimated SM background cross sections and the inoAMSBl benchmark 
point in fb for various multi-lepton plus jets +$ T topologies after cuts C2 with 
E C T = 100 GeV. 



and the inoAMSBl point, the total cross sections are shown in Table 4^ The hard 
C2 cuts prevent some backgrounds from producing noticeable cross sections indicated 
in the table by dashes. For a signal to be considered observable, we require that i.) 
S/B > 0.2, ii.) Signal has at least 5 events, and iii.) S > 5\/B (5cr significance), 
where S and B are respectively the signal and the background numbers of events for 
a given luminosity. 

Applying these criteria, the model's m.3/2 reach after 100 fb -1 of data collection is 
computed for each of the lepton channels and for each of = 100, 300 and 500 GeV. 



Figure 4.11 shows the plots of the cross section for the background (horizontal) and 
signal (decreasing slopes) as a function of m 3 / 2 . In the figure solid blue represents 
cuts with = 100 GeV, red dash-dotted represents 300 GeV, and purple dashed 
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Figure 4.11: Cross section for multi-jet plus $ T events with a). n(£) = 0, b). OS 
isolated dileptons c). isolated 3£s and d). isolated 4£s at the LHC after cuts C2 listed 
in the text with E C T = 100 GeV (blue solid), E C T = 300 GeV (red dot-dashed) and 
Ej, = 500 GeV (purple dashes), versus m 3 / 2 , from the inoAMSB model line points 
with tan/3 = 10 and fx > 0. We also list the 100 fkr 1 5a, 5 event, S > 0.1 BG limit 
with the horizontal lines. 



represents 500 GeV. Where the signal touches its corresponding (5a) background 
value on the plot is the where signal is no longer significant enough to be observable. 



These 777,3/2 values are listed in Table 4.4 (along with the SS dilepton reach). It is 
seen there that the best reach at 100 fb" 1 for the inoAMSB model (with tan (3 = 10 
and n > 0) is for multi-jet + $ T + 31 which shows significant cross section up to 
m 3/2 ~ 118 TeV (m- g ~ 2.6 TeV) when E C T = 500 GeV. 



HITs + multi-peak m(l + l ) Distribution 

We next examine the dilepton invariant mass distribution for cascade decay events 
for > 2 high p T jets, large $ T and SF/OS dileptons. It is known that the m(l + l~) 
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Table 4.4: Estimated reach of 100 fb" 1 LHC for m 3/2 (TeV) in the inoAMSB model 
line in various signal channels. 



distribution is capable of having a "kinematic mass edge" structure through the decay 
Z 2 -> or through Z 2 -> l ± l T Z 1 \gT\. Because the inoAMSBmodel predicts L/R 
splitting in the slepton masses, so it is expected that two edges can occur from the 
first type of decay through Z 2 — > lfl T and Z 2 — > 

It is also expected that a peak can occur in the m{l + l~) distribution around Mz 
because of the decays Z 3 — > ZZ\, Z 4 — ► ZZ\ and W 2 — >■ ZWi. In the case of point 
inoAMSBl, these branching occur for 25%, 6%, and 29% respectively. 



Figure 4.12 frame (a) shows the m(l + l~) distribution for inoAMSBl (red line) 
along with mAMSB (green) and HCAMSB(blue) points. The cuts are CI with $ T > 
300 GeV and A T > 900 GeV, and the SM backgrounds are completely suppressed. 
The peak at Mz is clear and there is a Z 2 decay double-edge structure calculated to 
be at values 



m~ / m% 

mtri-) < rn^/l - ^-y 1 - -4^ = 182 GeV and 304 GeV. (4.19) 



The former is for Z 2 to Ir decays where as the latter is for the decays to left-sleptons. 

Unlike inoAMSB, the mAMSBplot exhibits only a single edge because mAMS- 
Bpoints predict nearly-degenerate L and R sleptons. mAMSBalso lacks a Z-peak 
because the second heaviest neutralino, Z 2l is bino-like and thus has suppressed cou- 
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Figure 4.12: Dilepton invariant mass distribution, (a) comparison of inoAMSBto 
mAMSBand HCAMSB. CI cuts are used with $ T > 300 GeV and A T > 900GeV. 
(b) the same distribution except m 3 / 2 = 70 and 80 TeV values. The double-edge is 
still visible at these high m 3 / 2 values. 



plings to Z\Z . HCAMSBdoes not have an edge because here Z 2 is higgsino-like and 
has different decay products: Z 2 — > W^W T , Z-Ji, and Z\Z. Because the last decay 
channel should always be open, HCAMSBhas a smooth m(l + l~) distribution with a 
visible Z-resonance. 

We conclude that inoAMSB, mAMSB, and HCAMSBmodels can be distinguished 
at the LHC. First, AMSB models in general are identified by a wino-like Z\ with near 
degeneracy to the lightest chargino. This would be seen in gluino and squark cascade 
decay events including a long-lived chargino with a terminating track (HIT signal). 
Then the different m{l + l~) distributions could be used to distinguish between the 



inoAMSB, HCAMSB, or mAMSBcases. These results are summarized in Table 4.5 



for the reader's convenience. 
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Model 


HITs 


# of edges 


~%2 


? 

mr PS mr 


inoAMSB 


Yes 


2 


bino 


No 


mAMSB 


Yes 


1 


bino 


Yes 


HCAMSB 


Yes 





higgsino 


No 



Table 4.5: The three AMSB models can be distinguished in the dilepton invariant 
mass distribution by the number of edges. The structures depend crucially on the 
type of Z 2 and the slepton masses predicted by the models. 

4.5 Summary 

In this chapter, we have examined the phenomenology of supersymmetric models with 
the boundary conditions m ~ A ~ at M string , while gaugino masses assume the 
form as given in AMSB. We have labeled this class of boundary conditions gaugino- 
AMSB boundary conditions, or inoAMSBfor short. Such boundary conditions can 
arise in type IIB string models with flux compactifications. They are very compelling 
in that off-diagonal flavor violating and also CP violating terms are highly suppressed, 
as in the case of no-scale supergravity or gaugino-mediated SUSY breaking models. 
However, since gaugino masses assume the AMSB form at M GUT , the large U(l)y 
gaugino mass Mi pulls slepton masses to large enough values through renormaliza- 
tion group evolution that one avoids charged LSPs (as in NS or inoMSB model) or 
tachyonic sleptons (as in pure AMSB models). 

The expected sparticle mass spectrum is very distinctive. Like mAMSB and 
HCAMSB, it is expected that lightest neutralino (Zi) will be wino-like, and a quasi- 
stable chargino (Wi) would leave observable highly ionizing tracks in collider de- 
tectors. The spectrum is unlike mAMSB in that a large mass splitting is expected 
between left- and right- sleptons. The case when the string scale, M string , is much 



100 



4.5 Summary 



lower than Mqut was also investigated. In this case, the entire spectrum becomes 
somewhat spread-out, and if M s ~ 10 14 GeV, then the left-sneutrino becomes the 
LSP, which is excluded by double beta decay experiments [5]. 

We also saw in detail some aspects of LHC collider signatures. Since rriq < m- g 
in inoAMSB models, we expect dominant qq and qg production at LHC, followed by 
2-body q and g decays. This leads to collider events containing at least two very high 
Pt jets plus $ T as is indicative of squark pair production. 

While squark and gluino cascade decay events should be easily seen at the LHC 
(provided 7713/2 ~ 110 TeV), the signal events should all contain visible HITs, which 
would point to a model with ~ m % , as occurs in anomaly- mediation where 
M2 < Mi, M3 at the weak scale. We find a LHC reach, given 100 fb _1 of integrated 
luminosity, out to m 3 / 2 ~ 118 TeV, corresponding to a reach in m g of about 2.6 
TeV. We note here that if a signal is found at the outer edges of the reach limit, the 
signal will consist of typically just a few (5-10) events (due to hard cuts) over a small 
background. The signal events should include the characteristic presence of HITs, 
which should be absent in the background. As data accumulates, signals should also 
appear in the complementary channels, thus building confidence in a discovery. 

We also find that the invariant mass distribution of SF/OS dilepton pairs should 
have a distinctive two-bump structure that is indicative of neutralino decays through 
both left- and right- sleptons with a large slepton mass splitting. This distribution 
would help distinguish inoAMSB models from HCAMSB, where a continuum plus a 
Z-peak distribution is expected, or from mAMSB, where the two mass edges (present 
only if m is small enough that mj L and mj R are lighter than m^) would be very 
close together, and probably not resolvable. 
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Dark Matter in AMSB Models 



5.1 Introduction to Dark Matter in AMSB Models 

Now that we know what to expect for mAMSB, HCAMSB, and inoAMSBmodels at 
the LHC, we can now turn our attention to using astrophysical methods to observe 
them. In particular we are interested in i) whether these models can account for the 
observed cold dark matter (CDM) abundance 

VLh 2 = 0.1123 ±0.0035 68%CL (5.1) 

according to WMAP7 [75] measurements, where Q = and p c = 1.88 x 10~ 29 h 2 
g cm -3 is the critical closure density; ii) if current or future experiments can either 
confirm or rule out regions of their parameter spaces; and Hi) how cosmological data 
compares with LHC data. 

At first sight there appear to be problems with AMSB models arising from the 
fact that the LSPs annihilate and co-annihilate (with the lightest chargino) too effi- 
ciently [H] and leave the relic density of thermally produced neutralinos 1-2 orders 
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of magnitude below the experimental value in Equation 5^ Indeed, this is a major 
issue with AMSB models that has led many to consider other model types to be more 
interesting. However, it is argued here that when one looks back into the history of 
the Universe it is important not only to consider thermal production of LSPs in quasi- 
thermal equilibrium, but it is also important to consider the impact of non-thermal 
production of LSPs: production of LSPs through the decay of heavy particles that 
may have existed early on. 

In this chapter we will explore several possibilities to increase the CDM abun- 
dance^) that are strongly motivated in theory. String theories require the presence 
of heavy moduli fields (and other heavy scalars) and local SUSY theories require grav- 
itinos, both of which decay through gravitational interactions and can add to Qh 2 . 
We also consider here the case where the strong CP problem is solved by invoking the 
axion (a) solution. When this solution is combined with the SUSY framework, the 
axion's supersymmetric partner, the axino (a), appears in the spectrum [93] and its 
mass is tied (model-dependently) to the SUSY breaking sector. This means that the 
axino can be heavy enough that its decay chains terminate with LSP production or 
it can even serve as the LSP itself and contribute, along with the axion, to the CDM 
abundance. 

Non-thermal production mechanisms must not only produce enough DM to ac- 
count for WMAP, but must also be careful not to disrupt Big Bang Nucleosynthesis 
(BBN) which includes measurements of light element abundances, baryon-entropy 
ratio, and the neutron-proton ratio. With these requirements being met, the AMSB 
models' parameters can be further constrained and the reaches will be shown for di- 
rect and indirect detection experiments. 



This chapter is organized as follows. In Section 5.2 four non-thermal production 
mechanisms of dark matter will be described. The first two result in an increase of 
lowo-neutralino DM as a result of moduli decay and (thermally-produced) gravitino 
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decay. The remaining mechanisms will invoke the Peccei-Quinn (PQ) solution to the 
strong CP problem. The third mechanism will increase the wino DM through heavy 
axino decays and there will be an additional component of axion (a) DM from the 
breaking of the PQ symmetry. In the fourth mechanism we take the axino (a) to be 
the LSP such that it is produced both thermally and in neutralino decays. This last 
case will have a mixture of axion and axino CDM. 

Dark matter may be observable in future experiments through direct detection 
(DD) where the relic DM particles scatter from target nuclei in the experimental 
apparatus. It is also possible that it will be observed through indirect detection (ID) 
where DM annihilate into SM particles somewhere else in the galaxy and we see its 
products. A brief description of CDM experiments is given Section |5.3 



Sections 5.4| and 5.5 contain the results of this research. In Section 5.4, the cal- 



culated DD cross sections for AMSB wino DM are given. In Section 5.5, the ID 
rates for mAMSB, HCAMSB, and inoAMSBmodels are also given. The discussions 
in these sections describe the interactions that lead to the calculated rates and explore 
the current experimental exclusions and future experimental reaches of the parame- 
ter space. Finally, we close this chapter in Section [BTo" with a summary and conclusion. 



5.2 Thermal and Non-thermal Production of Dark 
Matter in AMSB Models 

As described in the Introduction of this chapter, dark matter (DM) production in 
the early universe can happen both thermally and non-thermally. In this section 
we explore four cases where the abundance of DM is composed of thermal and non- 
thermal components. In addition, we will see in these cases the DM can be wino (as 
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for our AMSB models), axion, axino, or mixtures of these. 

Neutralino CDM via Moduli Decay 

Shortly after the introduction of AMSB models, Moroi and Randall proposed a solu- 
tion to the AMSB dark matter problem based on augmented neutralino production 
via the decays of moduli fields in the early universe (EH]. The idea here is that string 
theory is replete with additional moduli fields: neutral scalar fields with gravitational 
couplings to matter. In generic supergravity theories, the moduli fields are expected 
to have masses comparable to m 3 / 2 . When the Hubble expansion rate becomes compa- 
rable to the moduli mass m^, then an effective potential will turn on, and the moduli 
field(s) will oscillate about their minima, producing massive excitations, which will 
then decay to all allowed modes: e.g. gauge boson pairs, higgs boson pairs, gravitino 
pairs, ■ ■ ■ . The neutralino production rate via moduli decay has been estimated in 
Ref. |85j . It is noted in Ref. [2] that the abundance- given by 

m-fc'~o.i x (^m (^X" (W) (m^i) 3 ' 2 (5 . 2) 

Zl VlOOGeV/V 9* J VHA m <t> J 

with cry = 3 x 10 -24 cm 3 /sec - yields nearly the measured dark matter abundance for 

wino-like neutralino annihilation cross sections and m<j> ~ 100 TeVQ These authors 

dub this the "non-thermal WIMP miracle" . 

A necessary condition for augmented neutralino production via scalar field decay 

is that the re-heat temperature of radiation Tr induced by moduli decays is bounded 

by Tr ~ 5 MeV (in order to sustain Big Bang Nucleosynthesis (BBN) as we know 

it), and Tr < Tf , where T/ is the freeze-out temperature for thermal neutralino 

production Tf Q ~ m^/20. If T R exceeds Tf Q , then the decay-produced neutralinos 

1 In inoAMSB models, we expect moduli with SUSY breaking scale masses, ~ rn^^/yV <C 
to 3 /2, where V is the (large) volume of the compactified manifold: V ~ 10 5 in Planck units. In this 
case, the mechanism would not so easily apply. 
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will thermalize, and the abundance will be given by the thermal calculation as usual. 

This "low re-heat" neutralino production mechanism has been investigated exten- 
sively by Gondolo, Gelmini, et. al. [58J[57J[56J. The low re-heat neutralino abundance 
calculation depends on the input value of Tr and the ratio b/m^, where b is the aver- 
age number of neutralinos produced in moduli decay, and is the scalar field mass. 
They note that theories with an under-abundance of thermally produced neutralino 
CDM with VL T ~ P ~ 10 -5 ( 100 GeV ) can always be brought into accord with the mea- 

^1 \ m z 1 ) 

sured DM abundance for at least one and sometimes two values of Tr^\ 

While the low Tr ~ 10 — 1000 MeV scenario with DM generation via scalar field 
decay is compelling, we note here that it is also consistent with some baryogenesis 
mechanisms: e.g. Affleck-Dine baryogenesis wherein a large baryon asymmetry is 
generated early on, only to be diluted to observable levels via moduli decay (69], or 
a scenario wherein the baryon asymmetry is actually generated by the moduli decay 

m 

Neutralino CDM via Gravitino Decay 

An alternative possibility for augmenting the production of wino-like neutralinos in 
AMSB models is via gravitino production and decay in the early universe. While 
gravitinos would not be in thermal equilibrium during or after re-heat, they still can 
be produced thermally via radiation off ordinary sparticle scattering reactions in the 
early universe. The relic density of thermally produced gravitinos as calculated in 
Ref's [12] [35] is given by 

- ±«* (i + 4) * ft) {Jkv) ■ M> 

2 References [56] [57] [58] also shows that an overabundance of thermally produced neutralino CDM 
can also be brought into accord with the measured abundance via dilution of the neutralino number 
density by entropy injection from the <f> field decay. Since this case does not attain in AMSB models 
(unless mg ~ 1300 GeV), we will neglect it here. 
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where gi and M { are the gauge couplings and gaugino masses evaluated at scale 
Q = T R , and u t = (0.018,0.044,0.117) and h = (1.266, 1.312, 1.271). Each gravitino 
ultimately cascade decays down to the wino-like Z\ state, so the neutralino relic 
density is given by 

_ 777,^ 

(5.4) 



h 2 = n T /h 2 + ^n T /h 2 . 



A plot of the value of Tr and m 3/ / 2 which is required to yield Q^h 2 = 0.11 from 



Equation 5.4 is shown in Figure 5.1 for mAMSB (mo = 0.01m 3 / 2 ), HCAMSB (a = 
0.02) and inoAMSB using tan /3 = 10 and \i > 0. The region above the f)g /i 2 = 0.11 
curves would yield too much dark matter, while the region below the curves yields 



too little. We should consider the curves shown in Figure 5.1 as only indicative of 
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Figure 5.1: Plot of allowed region of Tr vs. m 3 / 2 plane allowed for wino-like neutralino 
DM from thermal production plus thermally produced gravitino decay. For HCAMSB 
a = 0.02 and for mAMSB m = 0.01 m 3/2 . 



the simplest scenario for wino production via gravitino decay. Three other effects can 



substantially change the above picture from what is presented in Equation 5.4 



On the one hand, if moduli fields m exist with mass m^ m > 2m 3 / 2 , then 
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gravitinos can also be produced via moduli production and decay [73] [TU]. The 
exact abundance of these moduli-produced gravitinos is very model dependent, 
and depends on the moduli and gravitino mass and branching fractions. 

A second case arises if we consider gravitino production via inflaton decay at 
the end of inflation [HI]. This production mechanism depends on unknown 
properties of the inflaton: e.g. its mass and branching fractions, and the re- 
heat temperature generated by inflaton decay. These latter quantities are very 
model dependent. 

Additional entropy production generated via the inflaton, moduli and gravitino 



decays may also dilute the above relic abundance in Equation 5.4 



We will bear in mind that these possibilities permit much lower or much higher values 
of Tr and 7713/2 than those shown by the QzJ 1 " 1 — 0.1 contour of Figure 



5.1 



Neutralino CDM from Heavy Axino Decays Mixed with 
Axion CDM 

A third mechanism for increasing the wino-like relic abundance is presented in Ref. 
[43] . in the context of the PQMSSM. If we adopt the Peccei- Quinn (PQ) solution 
to the strong CP problem within the context of supersymmetric models, then it is 
appropriate to work with the PQ-augmented MSSM, which contains in addition to 
the usual MSSM states, the axion a, the i?-parity even saxion field s, and the spin— 
i?-parity odd axino a. The axino can serve as the lightest SUSY particle if it is lighter 
than the lightest i?-odd MSSM particle. The a and a have couplings to matter which 
are suppressed by the value of the PQ breaking scale f a , usually considered to be in 
the range 10 9 GeV ~ f a ~ 10 12 GeVjEQ. 

In Ref. [43J, it is assumed that m„ > , where Z\ is the LSP. In the AMSB 
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scenarios considered here, we will assume T R ~ 10 10 GeV, so as to avoid overproduc- 
tion of dark matter via gravitinos. With these low values of Tr, we are also below the 
axino decoupling temperature Ta-dcp = 10 11 GeV ( 1Q i2 a Gc y ) 2 (lr) ' so ^ ne ax i nos are 
never considered at thermal equilibrium However, axinos can still be produced 
thermally via radiation off usual MSSM scattering processes at high temperatures. 
The calculation of the thermally produced axino abundance, from the hard thermal 
loop approximation, yields [37] 

where g s is the strong coupling evaluated at Q = Tr and N is the model dependent 
color anomaly of the PQ symmetry, of order 1. Since these axinos are assumed quite 
heavy, they will decay to gg or modes, which further decay until the stable LSP 
state, assumed here to be the neutral wino, is reached. 

If the temperature of radiation due to axino decay (Tp) exceeds the neutralino 
freeze-out temperature Tf a , then the thermal wino abundance is unaffected by axino 
decay. If Tp < Tf a , then the axino decay will add to the neutralino abundance. 
However, this situation breaks up into two possibilities: i). a case wherein the axinos 
can dominate the energy density of the universe, wherein extra entropy production 
from heavy axino decay may dilute the thermal abundance of the wino-like LSPs, and 
ii). a case where they don't. In addition, if the yield of winos from axino decay is 
high enough, then additional annihilation of winos after axino decay may occur; this 
case is handled by explicit solution of the Boltzmann equation for the wino number 
density. Along with a component of wino-like neutralino CDM, there will of course 
be some component of vacuum mis-alignment produced axion CDM: thus, in this 
scenario, we expect a WIMP/axion mixture of CDM. 



109 



5.2 Thermal and Non-thermal Production of Dark Matter in AMSB 
Models 



Mixed Axion/Axino CDM 

In this case, we again consider the PQMSSM, as in the last subsection. But now, we 
consider a light axino with m„ < m Zv so that a is the stable LSP [2]. Here, the 
thermally produced wino-like neutralinos will decay via Z% — » aj, so we will obtain 
a very slight dark matter abundance from neutralino decay: Q^ TP = ^-f^/i 2 , 
since each thermally produced neutralino gives rise to one non-thermally produced 



(NTP) axino. We will also produce axinos thermally via Equation 5.5 Finally, we 
will also produce axion CDM via the vacuum mis- alignment mechanism[lj: Q a h 2 ~ 
4 ( io^gIv ) ®i ( we ^ a ^ e ^ ere ^ e i^it ial mis-alignment angle 9{ ~ 1). The 
entire CDM abundance is then the sum 

n a ~ a h 2 = n% TP h 2 + n T a p h 2 + n a h 2 . (5.6) 



In this case, the TP axinos constitute CDM as long as m„ ~ 0.1 MeV. The NTP axinos 
constitute warm DM for m a ~ 1 GeV [63] , but since their abundance is tiny, this fact 
is largely irrelevant. The entire CDM abundance then depends on the parameters f a , 
rria and Tr; it also depends extremely weakly on VL^Ji 2 , since this is usually small in 
AMSB models. 



As an example, we plot in Fig. 5.2 the three components of mixed axion/axino 
DM abundance from HCAMSB benchmark point 1 in Reference a = 0.025, 

7713/2 = 50 TeV, tan (3 = 10 and /i > 0. The neutralino thermal DM abundance would 
be ^zj 1 " 2 = 0.0015 if the Z\ was stable. We require instead fl a a,h 2 = 0.11, and plot 
the three components of VL a5 h 2 versus f a /N, for three values of Tr = 10 6 , 10 7 and 
10 8 GeV. The value of m„ is determined by the constraint fl a a,h 2 = 0.11. We see 
that at low values of f a /N, the NTP axino abundance is indeed tiny. Also the axion 
abundance is tiny since the assumed initial axion field strength is low. The TP axino 
abundance dominates. As f a /N increases, the axion abundance increases, taking 
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an ever greater share of the measured DM abundance. The TP axino abundance 
drops with increasing f a /N, since the effective axino coupling constant is decreasing. 
Around f a /N ~3 x 10 11 GeV, the axion abundance becomes dominant. It is in this 
range that ADMX[3] would stand a good chance of measuring an axion signal using 
their microwave cavity experiment. 



In Figure 5.3, we again require £l aa h = 0.11 for HCAMSB benchmark point 1, 




f/N [GeV] 



Figure 5.2: Abundance of TP and NTP axino DM and vacuum-misalignment pro- 
duction of axion CDM versus f a /N, for various values of Tr. 



but this time plot the value of Tr which is needed versus m„, for various values of 
fa/N. The plots terminate at high Tr in order to avoid reaching the axion decoupling 
temperature T a _d cp . Dashed curves indicate regions where over 50% of the DM is 
warm, instead of cold. Solid curves yield the bulk of DM as being cold. 

We see that for very light axino masses, and large values of /„, the value of Tr 
easily reaches beyond 10 6 GeV, while maintaining the bulk of dark matter as cold. 
Such high values of Tr are good enough to sustain baryogenesis via non-thermal 
leptogenesis[79]. although thermal leptogenesis requires Tr ~ 10 10 GeV[l0]. Since f a 
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is quite large, we would expect that the dominant portion of DM is composed of relic 
axions, rather than axinos; as such, detection of the relic axions may be possible at 



ADMX[3J. While Figures [53] and [573] were created for the HCAMSB model, quite 
similar results are obtained for the mAMSB or inoAMSB models. 




m. [GeV] 

Figure 5.3: Plot of Tr needed to ensure fl a ah 2 = 0.1 for HCAMSB benchmark Pt. 
1, versus m„ for various values of the PQ breaking scale f a . The dashed curves yield 
mainly warm axino DM, while solid curves yield mainly cold mixed axion/axino DM. 



5.3 Dark Matter Detection 

In this section a brief description of Dark Matter experiments are given. These break 
down into two categories: direct detection (DD) and indirect detection (ID) and each 
are described separately. 
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Direct Detection 

Direct detections experiments are those that are intended to measure properties of 
Weakly Interacting Massive Particles (WIMPs) via their direct elastic scatterings off 
of nuclei. This is in contrast to the ID experiments (described in the next subsection) 
where WIMP interactions lead to different final state particles, and these products 
are observed and analyzed. 

DD experiments are considered to be extremely important in WIMP searches for a 
number of reasons. Although a neutralino mass should in principle be re-constructable 
through LHC data, its interactions are much more difficult to decipher as a result 
of its inertness in the detector. If there are WIMPS with weak interactions, DD 
experiments should be able to probe the couplings of these particles. Furthermore, 
$ T signatures in a collider experiment do not give us information on the stability of 
particles on cosmological time scales, but only on collider time scales. Thus only DD 
experiments can tell us whether a particle produced in collisions is "absolutely" stable. 
And finally, DD experiments do not have the serious systematic uncertainties that ID 
experiments have. This is because the WIMP properties are measured directly, and 
because the backgrounds to these experiments are well understood and, in principle, 
can be controlled. The largest uncertainties in direct detection experiments comes 
from the observed local density of WIMPs in the vicinity of the earth and their 
velocity distributions. Increased local density of DM through clumping, etc. would 
have an impact on the rate of interactions in DD experiments. 

Dark matter must be cold and the calculations of wimp-nucleon scattering cross 
sections must be made in the zero momentum transfer limit. In SUSY theories, the 
interactions of neutralinos with quarks in the proton are dominated by t-channel CP 
even Higgs exchangers] . This is because the relevant vertex is gaugino-higgsino- 
Higgs [20]. Thus, purely gaugino or purely higgsino neutralinos will have suppressed 
3 s-channel squark exchange is also possible, but suppressed. 



113 



5.3 Dark Matter Detection 



couplings to nuclei [6 2J. 

The total elastic scattering cross section for WIMP-nucleon has contributions from 
a part that is spin- in dependent and from one that is spin-dependent. Because the 
WIMP de Broglie wavelength is expected to be of the order the dimension of heavy 
atoms, the WIMP will scatter from the coherent composite of nucleons. Thus the 
spin- independent cross section will be enhanced by factors of A 2 , where A is the 
atomic number of the target nucleus. The spin- dependent contribution to the cross 
section only occurs through spin interactions and scales with J(J + 1). 

The best limits on the spin-dependent cross section come from the Cryogenic Dark 
Matter Search (CDMS: Ge/Si detectors) experiment which exclude WIMPs with weak 
scale masses and cross sections at the level of 5 x 10 -8 pb level. We will compare our 
model to the spin-independent limits set by CDMS and other experiments in the text. 
We will also compare spin-dependent cross sections as well, but we should note that 
the constraints placed by spin- dependent DD experiments however are not nearly as 
strong as for spin- independent cross sections. 

The results of the XENON experiments (liquid Xenon detectors) are expected 
to surpass limits set by its competitors due to its size. XENON has already taken 
DM searches from the 10 to 100 kg-scale detectors. Impressively, future XENON 
experiments will be at the ton-scale with SI reaches in the 10~ 10 — 10~ n pb range. 

Indirect Detection 

Indirect detection (ID) observations are complementary to direct detection (DD) of 
DM, and there is a wide variety of ID programs that seek annihilation products 
of DM particles from displaced sources. Any locus of strong gravitation potential 
is capable of squeezing WIMPs close enough to interact highly with one another. 
Examples include the galactic center, dwarf spheroidal galaxies, the Sun, the core of 
the Earth, and inhomogeneities in WIMP halo profiles. However, different sources 
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will provide different information about DM depending on the strength and distance 
of the source, the distribution of DM around the source, and on the mass and coupling 
of the WIMPs. A few general descriptions of potential ID signals are described in 
this subsection^ 

Gamma Rays 

The galactic center of the Milky Way and dwarf spheroidal galaxies in or near the 
Milky Way are the most capable of gathering DM because of their large gravitation 
pull. But these sources are relatively far away and so evidence of WIMP interactions 
must arrive at Earth unattenuated and undeflected by large magnetic field variations 
in the galaxy. Because gamma rays do not carry charge and travel in straight lines 
they retain their spectral information and hence are excellent probes of interactions 
at far distances. 

The gamma ray fluxes produced by annihilations of DM for any structure and at 
any distance was performed in [32J and analyzed in [98]. The expression for the flux 
of gamma rays is 

- M * ir ' WV1 ^( 3xlO-^ m V 5 ) (^) V^)A«, (5.7) 

where AQ is the solid angle observed and the interaction cross section appears ex- 
plicitly and is the spectrum. J(Af2, ip) depends only on the DM distribution and 
is averaged over the solid angle to give 

^^(iG^)7/ (r(w)R <5 - 8) 

where the DM density, p, appears quadratically, ip is the angle from the galactic 
center, and the integral takes place over the line of sight (los). 
4 see [E] [62] for reviews. 
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The spectrum for standard astrophysical sources is generally expected to have 
power law behavior [34J, i.e. c ^r- oc E~ a . The HESS collaboration has found a high 
energy source of gamma rays near the galactic center jl] that exhibits the power law 
behavior with a = 2.2 ±0.09 ±0.05 above the threshold 165 GeV. It is difficult to pick 
out a signal of DM annihilations from this background, but it is still an interesting 
target for Fermi-LAT that explores E y > 1 GeV, which is well below the high HESS 
energy threshold. 

Antimatter 

Charged matter including positrons (e + ), antiprotons (p), and antideuterons (D) are 
also interesting probes of DM and are produced in WIMP pair-annihilations in the 
halo. In the case of wino DM pairs will annihilate most efficiently to combinations of 
Higgs and gauge bosons, whose further decays lead to antimatter [62J. Of these prod- 
ucts, cosmic positrons are the most interesting. Because positrons lose the majority 
their energy over kiloparsec lengths, they only probe the local DM (in sharp contrast 
to gamma rays). Since DM densities are better understood locally than, say, at the 
galactic center, positron measurements are subject to fewer systematic uncertainties 
than other antimatter species. 

Cosmic ray protons interacting with our galaxy's interstellar gas would provide 
the largest background sources of D and p fluxes. These types of reactions rarely pro- 
duce low-energy particles, much less nucleons with energies small enough to match 
the binding energy of D (2.2 MeV) jl9]. D production is indeed very rare and the 
background of these particles is highly suppressed at low energies, with estimated flux 
of $5 ~ 5 x 10 -13 [GeV cm 2 s sr] _1 . For the future GAPS[54J experiment, sensitive 
to fluxes above 1.5 x 10~ n [GeV cm 2 s sr] -1 , the background is essentially absent [43]. 
An anomalous amount of detected low-energy D could be a signature of neutralino 
annihilations in the halo, because the hadronization of WIMP hadronic annihilation 
channels {Z\Z\ — > qq) are confined by color to occur in the rest frame of the halo. 
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Vector boson and Higgs channels can also inject significant amounts of D for light 
Higgs and lower values of WIMP mass [15], especially near the D threshold. 

Antiproton production is not rare, in contrast to D, but it also would not be a 
clean signal of WIMP annihilation at low energies. This is because, of the two, p is 
more susceptible to ionization losses, synchrotron radiation, and solar modulation, 
allowing it to eventually populate the lower energies [49J. 

Neutrinos 

WIMPs are thought to become captured by elastically scattering from hydrogen and 
helium in the Sun and dissipating energy. High-energy neutrinos may be produced in 
WIMP annihilations in the core of the Sun. The z/^s can then undergo + q — > q' + fi 
scattering to produce muons in terrestrial ice or water. Thus, the detection of muons 
on Earth from the WIMP annihilation to neutrinos depends critically on the Sun 
capture rate of WIMPs, which further depends on the WIMP's nucleon cross section, 
its density near the Sun, and its mass. 

The the capture rate depends on the effective WIMP-nucleus cross section which 

is 

°*it = 4 D + <*% + 0-074 7 e , (5.9) 

where the factor 0.07 reflects the relative abundance if He to H along with other 
dynamical factors An important point here is that while spin-independent 

WIMP-nucleus cross sections tend to be highly constrained by DD experiments, SD 
cross sections are not so constrained and can be conceivable much higher. Thus ID 
experiments such as ICECUBE should be able to probe much lower values of the 
spin-dependent cross section than DD experiments such as COUPP. 
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5.4 Direct Detection Rates for the AMSB Models 

The ability to directly detect WIMPs depends not only on the mass, interactions, and 
the local density of the candidate, but also on the local velocity distribution. In the 
standard cosmological scenario WIMPs are thermally produced in an expanding uni- 
verse, they freeze out, and then fall into gravitational wells. It is generally assumed, 
at present, that WIMPs conform to a Maxwell-Boltzmann velocity distribution, i.e., 
f(v) ~ v 2 e~ v2 / v o, where ~ 220 km/s is the Sun's velocity about the galactic center. 
The same distribution holds for our cases, even despite the fact that non — thermal 
production of the DM candidates (moduli, gravitino, axino decays) can distort the 
initial velocity distributions of the WIMPs. This is because the original velocity 
distributions are red-shifted away, and the Maxwell-Boltzmann distributions arise as 
usual from gravitational in-fall. Thus, the calculations or WIMP-target scattering 
cross sections for direct detection can be carried out as usual. 



As described in Section 5^3 WIMP scattering cross sections from target nuclei are 
described through spin-independent (SI) and spin-dependent (SD) parts. Because SI 
cross sections receive enhancement from the mass of heavy targets DD will have bet- 
ter constraints of SI interactions. The SI WIMP-nucleon cross sections are calculated 
using the Isatools subroutine IsaReS [19] . 



In Figure 54 we scan over all parameters of the mAMSB, HCASMB, and inoAMS- 
Bparameter spaces, with a low and high value of tan /3 and \x > 0, and plot the Zi—p 
scattering cross section versus m^. Only points with > 91.9 GeV are retained, 
as required by LEP2 . The inoAMSBpoints appear as lines because they have one 
less parameter than either of the mAMSB or HCAMSB models. 

Several crucial features emerge from this plot. 

• First, we note that for a given value of , the value of a(Zip) is bounded from 
below, unlike the case of the mSUGRA model. That means that wino-WIMP 
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dark matter can be either detected or excluded for a given value. 

Second, we note that the cross section values generally fall in the range that 
is detectable at present or future DD experiments. The purple contour, for 
instance, exhibits the CDMS reach based on 2004-2009 data, and already ex- 
cludes some points, especially those at large tan (3. We also show the reach of 
Xenon-100, LUX, Xenon-100 upgrade, and Xenon 1 ton[53]. These experiments 
should be able to either discover or exclude AMSB models with values be- 
low ~ 90, 100, 200 and 500 GeV respectively. These WIMP masses correspond 
to values of m g ~ 690, 770, 1540 and 3850 GeV, respectively! The latter reach 
far exceeds the 100 fb _1 of integrated luminosity reach of LHC for trig, which 
were shown in Chapters [3] and [4] to be m g ~ 2.2 — 2.4 for HCAMSBand m g < 2.6 
TeV for inoAMSBrespectively. For inoAMSB models, where the minimal value 
of a SI (Zip) exceeds that of mAMSB or HCAMSB for a given mj> value, the 
Xenon 1 ton reach is to ~ 800 GeV, corresponding to a reach in m g of 6200 
GeV! 



In Figure [5T5] we scan the parameter spaces again and plot the spin-dependent DD 
cross section versus . The recent DD limit by the COUPP experiment appears 
at least two orders of magnitude above the theory predictions of all of the models. 
The probes of the SD cross section are actually better constrained by ID experiments. 
The ID limits from ICECUBE and its projected DeepCore limit are shown in this 
figure, and the cross section is that of winos interacting with solar Hydrogen. Though 
DeepCore will be able to access portions of parameter space, it will not be capable of 
reaching any of the the models' lower limits. 
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Figure 5.4: Spin-independent Z\ — p scattering cross section versus for mAMSB, 
HCAMSB and inoAMSB models for tan (3 = 10 and 40 and /x > 0. The parameters 
777.3/2 and also mo (for mAMSB) and a (for HCAMSB) have been scanned over. We 
also show the CDMS limit and projected Xenon and LUX sensitivities. 

5.5 Indirect Detection Rates for the AMSB Mod- 
els 

In this section the Indirect Detection (ID) rates for mAMSB, HCAMSB, and inoAMS- 
Bare given. We consider the indirect detection of WIMPs from their annihilation in 
the solar core via neutrino telescopes, and we also consider gamma ray and antimatter 
fluxes originating in halo WIMP annihilations. These fluxes depend (quadratically) 
on the assumed density profile of the galaxy, and the results of two profiles, isothermal 
and Navarro- Frenk- White [86] , are given. Most halo models are in near-agreement at 
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Figure 5.5: Spin-dependent Z\ — p scattering cross section versus for mAMSB, 
HCAMSB and inoAMSB models for tan (3 = 10 and 40 and /x > 0. The parameters 
/7i3/2 and also mo (for mAMSB) and a (for HCAMSB) have been scanned over. We 
also show the COUPP and IceCube limits in a SD (Zip). 



the Earth's position (~ 8 kpc from the galactic center), however they differ widely at 
the inner parsecs of the galaxy. This translates into large uncertainties in the gamma 



fluxes from these regions. As mentioned in Section |5.3[ antimatter signals tend to 
have fewer uncertainties because they should originate closer to earth. 

This section is broken into three parts. In the first subsection we consider the 
mAMSBmodel and calculate ID rates by varying itlq and tan and analyze the in- 
teractions that lead to these rates. In the following subsection we do the same for 
HCAMSB, but this time we vary a and tan (5. And finally, in the last subsection we 
show the calculated ID rates while varying the remaining parameter, 7713/2, for the 
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mAMSB, HCAMSB, and the inoAMSBmodels. 



Indirect Wino Detection for mAMSB 

We begin the analysis with mAMSB. We choose m 3 / 2 = 50 TeV and tan = 10, 
which gives a wino-like neutralino with mass — 144 GeV. All of the ID rates for 



these parameters, while varying m , are shown in Figure [576] and each is discussed in 
turn. 



In Figure 5.6 (a) we show the SI DD cross section for comparison, along with 
approximate reaches for XENON-10 and XENON-100. For higher m , the value of 
\n\ drops, much like it does in mSUGRA when the focus-point region is approached. 
Thus, when tuq is large, Z\ picks up a higgsino component that enhances the SI DD 
cross section since the relevant vertex is Higgs-higgsino-gaugino. For mo ~ 3500 GeV, 
it is seen that the Z\ becomes directly observable for XENON-100. 



Figure 5.6 (b) shows the first ID rate which is the flux of muons from the solar 
direction. To calculate the flux we used the Isajet/DarkSusy interface (see Chapter 
[T] for descriptions of these programs) and require > 50 GeV as required for ICE- 
CUBE. When ZiS annihilate in the core of the Sun they produce z/ M s that travel to 
Earth and convert in water/ice to muons via the charged-current interaction. Z\s are 
captured by the Sun in the first place through its spin-dependent interactions with 
Hydrogen that is mainly sensitive to Z* exchange. Thus, for larger values of m where 
the neutralino has more higgsino content, the flux is expected to be larger due to the 
Z-Zi-Z 1 coupling. 



In Figure 5.6 (c) we show the expected flux of gamma rays from the galactic core 
with E~f > 1 GeV as required by the Fermi Gamma-ray Space Telescope (FGST) 
experiment. The gamma ray fluxes are due to neutralino annihilations and depend 
mainly on the Z\Zi — > W + W~ annihilation cross section which occurs via chargino 
exchange. This is followed by W — > qq' — > 7r° — > 77 to produce the gammas. The 
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Figure 5.6: Wino CDM direct detection (a) and indirect detection (b)-(f) rates versus 
m in mAMSB for m 3 / 2 = 50 TeV, tan (3=10, and /i > 0. 



rates do not change with mo because ZiS remain wino-like and the chargino mass 
does not change significantly. As expected, the type of profile plays a significant role 
in gamma detection, and the predictions for the two differ by over an order of mag- 
nitude. Still, in both cases the fluxes exceed the approximate reach of FGST. 



Figures 5.6 (d)-(f) show antimatter fluxes for e + , p, and D respectively that come 
from halo annihilations. For e + and p there are observable rates predicted for Pamela 
[89] . and D observables rates might be seen by GAPS [54J. These elevated rates 
reflect the wino — wino annihilation into W + W~ cross section. 



Now the same calculations are done and plotted in Figure 5.7 with m 3 / 2 = 50 



TeV and p, > 0, but this time tan is increased 40. The SI DD cross section in frame 
(a) has the usual enhancement due to increased higgsino content of Z\ at larger m 
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values. Additionally, there is an enhancement for low m values as well. As described 
in Section 



5.3 



crjr depends on Higgs and squark exchange diagrams between Z\ 



and quarks and through loops via Higgs exchange between Z\ and gluons. Now that 
tan /3 is higher, the heavy Higgs mass is quite light, e.g., for m ~ 600 GeV the heavy 
Higgs mass is uih ~ 152 GeV. This results in a huge DD cross section at low mo, 
and regions with m ~ 900 GeV and tan (3 = 40 are already excluded by DD WIMP 
searches. 



Related to this, Figure 5.7 (b) shows an increase in the muon flux for low mo 
and tan (3 = 40. Usually the spin- dependent cross section is considered to be the 
most important part of Z\ capture at high tan in the sun because the SI cross 
section receives minimal enhancement from the target (H and He) masses. However, 
because of the lighter Higgs effects of the last paragraph at low m , the SI cross 
section plays a significant role in Z\ capture and hence muon detection on Earth. Of 
course, this region is already excluded by DD experiments. We also note the interest- 
ing "anti-resonance" effect in $ M for niQ ~ 850 GeV. Here, itla ~ 2m^, and on the 
Higgs resonance Z\Z\ — > bb is the dominant annihilation mode. The 6-decays modes 
produce a softer distribution of z/^s than from vector bosons leading to fewer muons 
passing the > 50 GeV requirement, and hence the sharp dip in flux. Finally, at 
large m we again see an increase in muon flux due to the increased higgsino compo- 
nent of Z\. 



Finally, gamma and antimatter fluxes versus mo are shown in Figure 5/7 in frames 
(c) - (f). The rates are again flat due to the constant Z\Z\ annihilation rate to vector 
bosons. The exception occurs at mo ~ 850 GeV where the A-resonance enhances 
the annihilation rate [18]. At high mo, e + and p fluxes drop slightly because Z\Z\ 
annihilations to bb become prominent relative to the usual vector boson products. 
Subsequent decays lead to softer e + and p, and hence less detectability above the 
E e + t p = 20 GeV threshold. In frame (f), the D rate actually increases with m since 
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Figure 5.7: Wino CDM direct detection (a) and indirect detection (b)-(f) rates versus 
m in mAMSB for m 3 / 2 = 50 TeV, tan /3=40, and fi > 0. 



the softer Ds produced are still in the low energy range of detection. 



Indirect Wino Detection for HCAMSB 

In this subsection results similar to those of the last subsection are shown for the 
HCAMSBmodel for m 3 / 2 = 50 TeV, but now we vary the a parameter. It was 
seen in Chapter [3] that a ~ corresponds to pure anomaly mediation, while larger 
a values give increasing Mi at the GUT scale. This affects the RGEs of sparticle 
masses through their hypercharge numbers, and this leads to a left-right split sparticle 
spectrum at the weak scale. This effect, combined with the effect of a large t-quark 
Yukawa coupling, leads to light, dominantly left top squark state t\ at large a. Raising 
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a also leads to a diminished \fi\ value, and this results in Z\ acquiring more of a 
higgsino component. 



Figure 5.8 (a) shows the SI cross section for direct detection. Reminiscent of mo 
in mAMSB, the cross section increases for larger values of a, and this is, of course, 
due to the larger higgsino component of Z\. 



In Figure 5.8 (b) is the muon flux versus a for the HCAMSBmodel. At low a the 
flux is quite low, but increases considerable at high a. The higgsino component at 
high a leads to unsuppressed couplings to Z. Since the Zi-nucleon SD cross section, 
is mainly sensitive to Z*-exchange, a higher flux of muons are seen on the Earth. 



In Figure 5.8, in frames (c)-(f) are the plots of the gamma ray flux and antimatter 
fluxes. The high rates are due to the large wino — wino — > VV annihilation cross 
sections and are relatively flat with a. At large a the annihilation to bb states is 
enhanced leading to smaller detection rates for e + and p (softer spectrum does not 
pass experimental threshold) and slightly increased D (by soft energy requirements). 



In Figure 5J3, we perform the same calculations as in Figure 5^ but now 
tan (3 is chosen to be 40. In frame (a), the spin- independent DD cross section is 
enhanced relative to the tan /3 — 10 case because of the now much lighter, heavy 
and pseudoscalar Higgs (H and A) and the increased 6-quark Yukawa coupling. As a 
increases, m# increases and the rate diminishes until the highest a values are reached. 
Here the presence of the higgsino component of Z\ enhances the cross section. As in 
the mAMSBcase, the large SI enhancement leads to greater solar capture rate of Zis, 
which in turn leads to greater muon fluxes on earth, and this is seen in frame (b). 



In Figure 5.9 (b), we see the same type of anti-resonance effect as in Figure 5.7 



where near 2m % ~ tua, annihilation to heavy fermion pairs dominate over the usual 
VV pairs (V = W , Z), and lead to a softer distribution of neutrinos. In frames 
(c)-(f) are the halo annihilation rates for HCAMSBwith tan (3 = 40. These rates are 
generally flat with changing a, and do not suffer an increase as compared with low 
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Figure 5.8: Wino CDM direct detection (a) and indirect detection (b)-(f) rates versus 
a in H CAM SB for m 3/2 = 50 TeV, tan (3=0, and /x > 0. 



tan (3, since wino — wino — > VV still dominates the annihilation rate. The exception 
again occurs near the A-resonance, and halo annihilation is enhanced by pseudoscalar 
Higgs exchange. 



Indirect Detection Rates for mAMSB, HCAMSB, and inoAMS- 
B versus m 3 / 2 

In this subsection we explore the remaining parameter of our AMSB models: 7713/2. 
While m 3 / 2 is roughly bounded from below by LEP2 constraints, in principle there is 
no upper bound for the parameter. In general, increasing the gravitino mass increases 



all sparticle mass as can be seen for example in Table 4^ for the inoAMSBcase. Thus 
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Figure 5.9: Wino CDM direct detection (a) and indirect detection (b)-(f) rates versus 
a in HCAMSBfor m 3/2 = 50 TeV, tan 0=40, and /i > 0. 



we expect a diminution of rates toward higher m 3 / 2 values due to decreased interac- 
tions. Also, in all ID plots we assume the standard NFW halo profile. 



In Figure 5.10 we show direct and indirect wino DM detection rates versus 7B3/2 



for all three models considered in this thesis, with tan /5 = 10 and /i > 0. For the 
comparison, in the HCAMSBcase the choice a = 0.1 is made, and for mAMSBwe 
take m = 1 TeV. 



The spin- independent cross section is shown in frame (a) of Figure 5.10 for all 



models. The larger inoAMSBcross section is due in part to the smaller /x values for 
a given 777.3/2, which enhances Z\ scattering via Higgs exchange made possible by the 
Higgs-higgsino-gaugino vertex. 



Figure 5.10 (b) shows the muon flux with the minimum energy of > 50 GeV re- 
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quired for detectability at IceCube. Due to the relative values of /J in all three models, 
inoAMSByields the highest rates while mAMSByields the lowest. Since inoAMSB- 
generally has the lowest |/i| values, its Z\ can interact more strongly through Z* 
exchange (see previous two subsections). This enhances the SD cross section, which 
eventually leads to greater /x flux on earth. The rough reach of IceCube is shown, 
and the figure indicates that the lower values of m.3/2 may be accessible in — > \i 
searches. 



The ID rates for 7,e + ,p, and D are shown in Figure 5.10, frames (c)-(f). The 
rates for all models are nearly identical in each case for varying 7723/2. This is due to 
the dominance of Z\Z\-^- VV halo annihilations, which mainly depend on the gaug- 
ino (wino) component of Z\. Rough reaches for Fermi-LAT, Pamela, and GAPS are 
shown for reference, and high wino annihilations should yield observable signals. 



Increasing tan (3 to 40, the same rates are shown versus 7713/2 in Figure 5.11 



As in the previous subsections, the SI cross section in frame (a) is enhanced relative 
to tan (3 = 10 due to the simultaneous decrease in the Higgs mass and increase in 
6-quark Yukawa coupling. Because inoAMSBhas both the smallest m# and smallest 
|/i I values of all models for low m 3 / 2 , this model shows the largest SI cross section. 
As m.3/2 increases, m# increases for inoAMSBand HCAMSB, but actually decreases 
for mAMSB. Thus for 7773/2 ~ 75 TeV, the mAMSBmodel yields the highest value of 
o sl {Z lV ). 



Figure |5.11 frame (b), shows the muon flux, which is elevated relative to the 
tan (3 = 10 case, again due to the increased SI cross section of frame (a). At low 777.3/2 
the inoAMSBmodel shows the highest flux due to its low fi value which enhances the 
Z* exchange in Z\-q scattering. Around 7723/2 ~ 60 TeV, the mass of A is near twice 
the Z\ mass, and resonant Z\ annihilations to bb final states dominate over VV. This 
diminishes the flux of detectable muons (see previous two subsections for further dis- 
cussion). For higher 777.3/2 values, beyond the resonance, the mAMSB model yields 
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tan(3 = 10 




Figure 5.10: DD and ID rates for wino CDM in mAMSB, HCAMSB, and inoAMSB- 
models versus m 3/ / 2 , for tan (3 = 10, and \x > 0. For mAMSB, m = 1 TeV, and for 
HCAMSB, a = 0.1. In these plots we adopt the NFW DM halo profile. 



the highest muon flux due to its elevated value of a (Zip) 



Figure 5.11, frames (c)-(f), show the gamma and antimatter fluxes versus 7723/2. 
Relatively little change in the ID rates are seen in going from tan /? = 10 to tan = 40, 
since wino — wino — > VV dominates the annihilations. An exception, of course, is 
due to A-resonance enhancement of halo annihilations, where 2m % ~ tua,- And D 
rates appear enhanced for inoAMSBand mAMSBat large because in these cases 
the ZiZi—>bb annihilation rate, which does receive tan /3 enhancement, contributes 
to the detection of rather low energy Ds. 



130 



5.6 Summary 



tanfl = 40 



10 

10" 

l(f 

10"' 



1 


M 1 I | I 1 


i i i i i i i 

mAMSB 1 




\ 


HCAMSB^ 






■ — inoAMSB - 










(a) 




"l 


i i 


1 1 1 1 [ 1 [ 1 1 1 1 1 1 1 1 1 ] 




10" 

io- u 



10"' 



10" 



E > 1 GeV (Ferrai-LAT) 



10" 




80 100 120 140 160 





10 


'e 






L0 1 


□ 




y 






10° 



-I J I I I I I I I I I I I I I I I I I I I I I I I- 
S s (b) l 

*E > 50 Ge VffesCufe^ 

i I I I I I I I I I I I 




(■ 1 




z 1 1 11 


i 1 i i i | i i i | i i i | 


mm 


B 

u 
> 

o 


io" 7 




J 


(d) ^ 


t- 

© 


g 
— 


io s 


-£ + 

— f 


= 20 GeV (F^Ml 


e^A) - 


3 




-i 1 ii 


i i i i i 1 i i 


i i i i i r 



™ m (TeV) 




40 60 80 100 120 140 160 
^(TeV) 



Figure 5.11: DD and ID rates for wino CDM in mAMSB, HCAMSB, and inoAMSB- 
models versus m 3/ / 2 , for tan j3 = 40, and /i > 0. For mAMSB, m = 1 TeV, and for 
HCAMSB, a = 0.1. In these plots we adopt the NFW DM halo profile. 



5.6 Summary 

In this Chapter we have investigated aspects of cold dark matter in three models of 
anomaly mediation: mAMSB, HCAMSB and inoAMSB. Typically, each gives rise to 
a wino-like lightest neutralino, unless very high values of m (for mAMSB) or a (for 
HCAMSB) are used, in which case the Z\ becomes a mixed wino-higgsino state. In 
this class of models with a wino-like Z\, the thermal abundance of neutralino CDM 
is well below measured values, unless rug ~ 1300 GeV. We discuss four ways to 
reconcile the predicted abundance of CDM with experiment: 

1. enhanced neutralino production via scalar field (e.g. moduli) decay, 
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2. enhanced neutralino production via gravitino decay, where gravitinos may arise 
thermally, or by moduli or inflaton decay, 

3. enhanced neutralino production via heavy axino decay, and 

4. neutralino decay to axinos, where the bulk of CDM comes from a mixture of 
vacuum mis-alignment produced axions and thermally produced axinos. 

Cases 1 and 2 should lead to a situation where all of CDM is comprised of wino- 
like WIMPs; they will be very hard, perhaps impossible, to tell apart. Case 3 would 
contain a mixture of axion and wino-like WIMP CDM. It is a scenario where it is 
possible that both a WIMP and an axion could be detected. Case 4 predicts no 
direct or indirect detection of WIMPs, but a possible detection of relic axions. It is 
important to note that more than one of these mechanisms may occur at once: for 
instance, we may gain additional neutralino production in the early universe from 
moduli, gravitino and axino decay all together. 



In Sections |5 .4| and [5. 5[ we presented rates for direct and indirect detection of relic 
wino-like WIMPs. The SI direct detection cross sections are bounded from below. 
Ultimately, ton-scale noble liquid or SuperCDMS experiments should probe out to 
~ 500 GeV, which would exceed the 100 fb" 1 reach of LHC; a non-observation 
of signal would put enormous stress on AMSB-like models as new physics. We also 
evaluated SD direct detection: current experiments have little reach for AMSB-like 
models, although IceCube DeepCore and possibly COUPP upgrades may probe more 
deeply. 

WIMP indirect detection rates for all three AMSB models were also presented. 
The IceCube experiment has some reach for WIMPs from AMSB models, especially 
at high tan (5 or when the Z\ picks up a higgsino component. We noted an interesting 
inverse resonance effect in the muon flux detection rate, caused by transition from 
solar core annihilations to VV states, to annihilations to mainly bb states. The 
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detection of 7s, e + s, ps and Ds are all elevated in AMSB-like models compared to 
mSUGRA |18j . due to the high rate for Z\Z\ — > VV annihilation in the galactic halo. 
The results do depend on the assumed halo profile, especially for 7 ray detection in 
the direction of the galactic core. Generally, if a signal is seen in the e + channel, 
then one ought to be seen in the p channel, and ultimately in the 7, D (if/when 
GAPS flies) or direct detection channel. In addition, a sparticle production signal 
should ultimately be seen at the LHC, at least for rrig ~ 2400 GeV, once 100 fb _1 of 
integrated luminosity is accrued. 

As a final remark, we should note here that the dark matter detection signals 
all provide complementary information to that which will be provided by the CERN 
LHC. At the LHC, each model- mAMSB, HCAMSB and inoAMSB- will provide a 
rich assortment of gluino and squark cascade decay signals which will include multi- 
jet plus multi-lepton plus missing Et events. In all cases, the wino-like lightest 
neutralino state will be signaled by the well-known presence of highly ionizing tracks 
(HITs) from quasi-stable charginos with track length of order centimeters, before 
they decay to soft pions plus a, Z%. It was discussed in Chapter [i] that the three 
models should be distinguishable at the LHC by the differing opposite-sign/same 
flavor dilepton invariant mass distributions. In the case of mAMSB, with mi ~ m ; ~ , 
we expect a single mass edge from Z 2 —> Ul,r, Zi decay. In HCAMSB, the 

sleptons are rather heavy, and instead Z 2 — > Z\Z occurs at a large rate, leading 
to a bump in m{l + l~) ~ Mz, upon a continuum distribution. In inoAMSB, with 
m^ 2 > m i L R i but with II and Ir split in mass (due to different U(1)y quantum 
numbers), a characteristic double mass edge is expected in the m{l + l~) invariant 
mass distribution. 
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Conclusion 



We conclude by summarizing the main points of this thesis. We have compared three 
model types, all having supergravity anomaly contributions (AMSB) to the soft pa- 
rameters. We began with the mAMSBmodel of Chapter [2] which was the starting 
point of comparison for all other models. This model features a wino-like LSP (Zi), a 
lightest chargino iWi) with mass nearly equal to that of the LSP, and a near left-right 
degeneracy in scalar masses. Since ~ m wv we ex P ec t the presence of HITs (for 
all of the models) as described in Chapter [2J The second lightest neutralino, Z 2 , is 
bino-like and can be produced in LHC pp collisions and subsequently decays through 
Z2 — > li/^-i while slepton decays produce the LSP and dilepton final states. Since 
left- and right-handed states are degenerate, the OS m(/ ± / =F ) distribution for these 
final states results in a single (comparatively large) mass edge, as can be see in Figure 

We also examined the Hypercharged Anomaly Mediated Supersymmetry-Breaking 
model, which incorporates effects from a geometrically separated hidden sector. The 
hidden sector communicates the supersymmetry-breaking to the visible sector via a 
combination of supergravity anomaly and a new U(l) mediation. The soft masses are 
those of mAMSB, but with an extra contribution to the bino mass parameterized as 
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a in Equation (3.12). Because left- and right-handed particles have different hyper- 
charge quantum numbers, the mass spectrum for this model is highly left-right split 
(in contrast to mAMSB). The LSP is still a wino-like neutralino with near-degeneracy 
to the lightest chargino, and we again expect highly-ionizing track chargino events. 
The second-lightest neutralino, however, is higgsino-like in this case due to heaviness 
of the bino soft mass, M\. The higgsino component in Z 2 allows for decays to Z and 
Z\. Without intermediate slepton decays, the OS dilepton distribution appears as a 
smooth distribution with a peak around mz- 

The last model (s) we examined originate in string theories with Calabi-Yao orien- 
tifold compactifications and moduli stabilization through fluxes and non-perturbative 
(KKLT) effects. The supersymmetry-breaking is communicated to the MSSM via 
gravitation. The resulting soft contributions appear as a combination of anomaly 
mediation and gaugino mediation, and we have labeled these parameters space points 
"inoAMSB". In these models, the AMSB contributions are only present for gaugino 
masses, Mi,M 2 , and M3. But like gaugino mediation, the soft scalar and trilinear 
coupling are suppressed, i.e., m ~ A ~ at M str i ng . We expect for this class of 
models to yield HITs, as usual for AMSB models. We also expect with a double 
edge structure in the m(Z ± Z =F ) distribution. This is because, like in mAMSB, Z2 is 
bino-like, but the sleptons here are left-right split leading to two distinct mass edges, 



instead of just one as in the mAMSB framework (again, see Figure 4.12). 

Finally, we explored the Dark Matter predictions for all of these models. It is 
well-known that the neutralino LSPs of AMSB models annihilate and co-annihilate 
too efficiently to account for DM relic density of the universe. However, if any of 
these models are the correct description of nature, there are likely other non-thermal 
sources of LSP-production that can serve to increase the value of the relic density. 
Assuming these mechanism provide the correct DM abundance, we calculated direct 
and various indirect detection rates for a multitude of experiments for every model. 
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For the direct detection, we saw that CDMS data has already excluded some regions 
of the parameters spaces of these models. It was also striking to find that the cross 



sections were all bounded from below (Figure 5.4), leading to a XENON 1 Ton reach 
that far exceeded the reach of the LHC; the extreme example being the inoAMS- 
Breach nig ~ 6.2 TeV compared to m g ~ 2.6 TeV for 100 fb _1 of LHC data. 

For indirect detection, we found that the ICECUBE experiment provided some 
reach for AMSB models, particularly in cases with high tan /3, or in parameter space 
regions where Z\ has significant higgsino content (high mo or high a). The detection 
elevated rates of gammas, positrons, antiprotons, and antideuterons in AMSB models 
is due to high rate of Z\Z\ — > VV annihilations. We noted interesting inverse reso- 
nance effects in the muon flux. We also observed how indirect detection rates depend 
on the assumed galactic DM profile. 
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